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Abstract. We have analyzed deep Near Infrared (NIR) broad band images for a sample of Blue Compact Dwarf Galaxies 
(BCDs), observed with the ESO NTT* and Calar Alto** 3.6m telescopes. The data presented here allows for the detection and 
quantitative study of the extended stellar low-surface brightness (LSB) host galaxy in all sample BCDs. NIR surface brightness 
profiles (SBPs) of the LSB host galaxies agree at large galactocentric radii with those from optical studies, showing also 
an exponential intensity decrease and compatible scale lengths. At small to intermediate radii (within 1-3 exponential scale 
lengths), however, the NIR data reveals for more than one half of our sample BCDs evidence for a significant flattening of the 
exponential profile of the LSB component. Such profiles (type V SBPs, Binggeli & Cameron 1991) have rarely been detected in 
the LSB component of BCDs at optical wavelengths, where the relative flux contribution of the starburst, being stronger than in 
the NIR, can readily hide a possible central intensity depression in the underlying LSB host. The structural properties, frequency 
and physical origin of type V LSB profiles in BCDs and dwarf galaxies in general have not yet been subject to systematic studies. 
Nevertheless, the occurrence of such profiles in an appreciable fraction of BCDs would impose important new observational 
constraints to the radial mass distribution of the stellar LSB component, as well as to the photometric fading of these systems 
after the termination of star-forming activities. We test the suitability of two empirical fitting functions, a modified exponential 
distribution (Papaderos et al. 1996a) and the Sersic law, for the systematization of the structural properties of BCD host galaxies 
which show a type V intensity distribution. Either function has been found to satisfactorily fit a type V distribution. However, it 
is argued that the practical applicability of Sersic fits to the LSB emission of BCDs is limited by the extreme sensitivity of the 
achieved solutions to, e.g., small uncertainties in the sky subtraction and SBP derivation. We find that most of the sample BCDs 
show in their stellar LSB host galaxy optical-NIR colors indicative of an evolved stellar population with subsolar metallicity. 
Unsharp-masked NIR maps reveal numerous morphological details and indicate in some cases, in combination with optical 
data, appreciable non-uniform dust absorption on a spatial scale as large as ~1 kpc. 
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1. Introduction 

The star-formation history and chemodynamic evolution of 
Blue Compact Dwarf (BCD) galaxies are central issues in the 
contemporary dwarf galaxy research. In spite of being old 
in their vast majority, BCDs resemble in many aspects un- 
evolved low-mass galaxies in the early Universe. They are 
gas-rich (Hi mass fraction of typically >30%) and metal- 
deficient (Zq/SO^ Z <Zq/3) extragalactic systems, undergo- 
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ing intense, spatially extended star-forming (SF) activity. Such 
properties are believed to have been common among young 
low-mass objects at high to intermediate redshift, such as pre- 
galactic building blocks (Lowenthal et al. 1997 Hirashita et 
al. 120001 Fujita et al. 1200 1> or the progenitors of the present- 
day dwarf spheroidals (e.g. Babul & Rees 1992' Guzman et 
al. I1998> . BCDs are therefore convenient nearby laboratories 
to study at high spatial resolution the impact of collective star 
formation on the spectrophotometric and chemodynamic prop- 
erties of these distant and faint extragalactic sources. Moreover, 
they are important testbeds for deducing constraints to cosmo- 
logical parameters, such as the primordial ''He abundance ratio, 
and to monitor the synthesis and dispersal of heavy elements 
in a nearly pristine environment (Peimbert & Torres-Peimbert 
[T974I Pagel et aL [T992l Izotov et aL a997> . 

The understanding of the origin and implications of the 
starburst phenomenon in BCDs is necessary for elucidating 
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evolutionary pathways of dwarf galaxies (DGs) in general. Are 
BCDs active phases in the lifetime of dormant dwarf iiTeg- 
ulars (dis) and do the latter fade to dwarf ellipticals (dEs) 
once their gas reservoir has been depleted (see e.g. Lin & 
Faber [T983l Thuan [T985l Silk et al.[T987 Davies & Phillipps 
fl988l? What is the role of the environment (e.g.. Babul & 
Rees .1992. Pustil'nik et al. 12001^ and of Dark Matter (DM; 
Dekel & Silk [T986I Ferrara & Tolstoy l2000l . and does the 
latter invariably dominate the mass within the Holmberg ra- 
dius of a BCD (Papaderos et al. 1996b; hereafter P96bi) ? Do 
most BCDs undergo intermittent bursts or rather prolonged pe- 
riods of elevated star formation (Vallenari & Bomans 119961 
Noguchi .2001 Rieschick & Hensler 2001 Schulte-Ladbeck et 
al. l2001> ? Despite much previous effort, the observational ev- 
idence available thus far is still too fragmentary to allow for 
unambiguous answers to the aforementioned questions. 

Recent studies suggest, however, that key information for 
assessing DG evolution can be inferred from studies of the stel- 
lar low-surface brightness (LSB) host galaxy of these systems. 
In BCDs, the LSB component, underlying the SF regions, has 
first been disclosed through deep CCD imaging by Loose & 
Thuan (1986, hereafter LT86l, and has in the following been 
confirmed and further studied by various authors (cf., e.g., the 
list given later in this Section). This extended stellar host was 
found to account for ^1/2 of the light inside the 25 B mag/D" 
isophote (P9^ Salzer & Norton 1993, ™d to typically domi- 
nate the intensity and color distribution of BCDs for fi >24.5 B 
mag/n". Such an evenly distributed, evolved stellar component 
is observed in all types of DGs, except for the extremely rare 
type of iO BCDs in the classification scheme of LT86' Its red 
colors and smooth appearance in the main class of lE/nE BCDs 
( ILT86> indicate that these systems are several Gyr old, gas-rich 
DGs, having not been forming stars at the presently large rate 
throughout their lifetime. Different lines of evidence, outlined 
in the following, suggest that elaborate studies of the structural 
and kinematic properties of the LSB component are fundamen- 
tal to assess at least two central issues of DG research: the evo- 
lutionary connections between DGs and the regulation of the 
SF process in these systems. 

According to the standard evolutionary hypothesis for 
dwarf galaxies, dIs, dEs and BCDs differ basically by their gas 
content and the amplitude of their ongoing SF activity (Thuan 
Il985. Davies & Phillipps 1988 1. One would therefore expect 
that, on average, the evolved stellar LSB host in all these three 
main DG classes is indistinguishable from one another with 
respect to its structural properties. However. .P96b. and subse- 
quent authors (Patterson & Thuan 1 19961 Marlowe et al. 119971 
Salzer & Norton 1999 see also Papaderos et al. "2002 1 found 
that, at equal absolute magnitude, the stellar LSB component of 
iE/nE BCDs is systematically more compact than other types 
of DGs. P96b' interpreted this structural disparity as the result 
of adiabatic contraction of the stellar LSB component of BCDs 
in response to a large-scale gas inflow, preceding the ignition 
of a starburst in a BCD. Quite interestingly, subsequent in- 
terferometric H I studies have shown that the gaseous halo of 
BCDs is by at least a factor of two more centrally concentrated 
than in dIs (van Zee et al. 1998"200r, see also Simpson et al. 
[2000j . lending circumstantial support to the latter hypothesis. 



Whether or not such observational constraints can be accounted 
for by the dynamical evolution of the stellar and gaseous matter 
in DGs (lP96bt or, alternatively, an extraordinarily dense Dark 
Matter halo, being particular to BCD galaxies (Meurer et al. 
1998 1 awaits to be investigated by numerical 3D-simulations. 
These are also needed to explore a possible connection between 
the LSB morphology and the evolutionary state of BCDs, as 
proposed by Noeske et al. (I2000> . 

The evolved LSB component contains, due to its high 
M/L, the bulk of the stellar matter in a BCD. Thus, provided 
that DM does not dominate within the optical radius, it forms, 
together with the gaseous component, the graviational potential 
within which SF activity occurs. It is therefore worth explor- 
ing whether the structural properties of the LSB component in- 
fluence the SF process in a BCD. P96b found a trend for the 
fractional surface area of the SF component of BCDs to de- 
crease with increasing LSB luminosity. For a constant M/L, 
this trend translates into a mass-morphology connection for 
BCDs: SF activity in more massive BCDs occurs mainly in the 
inner part of the LSB component, leading to a nE morphology 
and surface brightness profiles (SBPs) sometimes superficially 
resembling a de Vaucouleurs law. Conversely, SF activity in 
low-mass BCDs is spread over a larger portion of the LSB host, 
resulting in an iE morphology and SBPs possessing a conspicu- 
ous plateau feature (cf. Papaderos et al. 1996a, hei'eafter .P96aj 
at intermediate intensity levels. 

Our current knowledge of the nature of the underlying LSB 
component relies mainly on optical surface photometry stud- 
ies of its fa int per iphery (e.g. lLT86l Ku nfli et aL [T9881 |P96al 
Telles et al. TW Marlowe et al. ^9971 DoubUer et al. [T9971 
.1999. Salzer & Norton 1999. Vennik et al.|2000, Cairos et al. 
l2001all200Tbl Makarova et al. l2002l . In optical wavelengths, 
extended stellar and gaseous starburst emission overshines the 
LSB component out to a galactocentric radius of typically ~2 
exponential scale lengths. Only at larger radii, where the star- 
burst emission is in most cases negligible, the LSB host can be 
studied directly. However, to explore its possible dynamical in- 
fluence on the global SF process, it is essential to pin down its 
intensity distribution at smaller radii, if possible just beneath 
the SF regions. Deprojection of SBPs would then allow one to 
put constraints on the density profile and the gravitational po- 
tential of the evolved stellar host (cf., e.g., P96ai. In addition, 
one would be able to correct optical colors inside the SF com- 
ponent for the line-of-sight contribution of the underlying stel- 
lar background (cf. Cairos et al. I2002al Papaderos et al. l2002> . 
thus better constrain the SF history of these systems. 

One way to alleviate the problem of the extended starburst 
emission is to extend studies to the Near Infrared (NIR). At 
these wavelengths, the young stellar populations and the ion- 
ized gas contribute a smaller fraction of the total light of the 
galaxy than in the optical. For instance, evolutionary synthe- 
sis models by Kriiger et al. ( 1995 ) predict that a moderately 
strong burst during its peak luminosity accounts for only ^^20% 
of a BCD's emission in the K, but for >80% in the B band. 
Observations of BCDs also show that starburst emission is 
weaker in the NIR, and becomes negligible at a smaller galacto- 
centric distance than in the optical (e.g. Vanzi et al. 1996" 2 0021 
Beck et al. .l997. Alton et al. .l994. James .l994j . NIR data al- 
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low therefore to study the stellar LSB component at smaller 
radii and, using optical-NIR colors (e.g. B — J), better con- 
strain its formation history. To achieve these objectives one 
needs to unambiguously detect the LSB component and study 
its intensity over a sufficient span (A/i > 2 mag). Empirical 
estimates, based on published optical and NIR data, suggest 
that these requirements are met by extending NIR studies to 
^22-24 J mag/n". So far, only a few BCDs have been studied 
at those surface brightness levels with an accuracy high enough 
to pinpoint structural properties and colors (Mkn 86, Gil de Paz 
et al. l2000al Tol 0645-376, Doublier et al. 200 ij Tol3, Vanzi et 
al.|2002|l. 

The present analysis is the first part of an observational 
project (see also Cairos et al. 2003; hereafter C03 ), aiming at 
a systematic study of the NIR properties of nearby BCDs with 
large array detectors on 4m-class telescopes. For this purpose, 
we take advantage of a large set of imaging data, homogeneous 
with respect to its limiting surface brightness and the methods 
used in its processing. The observations have been conducted 
so as to permit NIR surface photometry out to a comparable 
radius as in optical wavelengths, with the purpose of a multi- 
wavelength investigation of the LSB component. 

This paper is structured as follows: In Sect. 2, we describe 
the sample selection and data acquisition, and discuss the data 
reduction, photometric transformations and extinction correc- 
tions. Section 3 focusses on the derivation of SBPs and their 
decomposition into the luminosity components owing to the 
old LSB host and the younger stellar populations and SF re- 
gions. In Sect. 4, individual objects are presented in detail. Our 
results are discussed in Sect. 5. Sect. 6 summarizes this work. 



2. Observations and data reduction 

2. 1. Sample selection 

Our sample covers the whole morphological spectrum of 
BCDs, containing both examples of the main morphological 
class (iE/nE systems according to ILT86t . and of the less fre- 
quent iI/iI,C/iO BCD classes. The latter subset includes the 
metal-poor galaxies Tol 65 and Tol 1214-277 (Z <Zq/20). 
The sample also includes one intrinsically luminous blue com- 
pact galaxy (BCG), UM448. 

Table [Olists the adopted distance to each BCD. Whenever 
available, literature distances based on standard candles were 
preferred for nearby objects. Those literature distances that 
rely on redshifts have been corrected for peculiar motions 
within the local supercluster and assume a Hubble constant 
Hq =75 km s^^ Mpc~^. When no literature data was available, 
the distance was inferred from the Hq and the heliocentric ve- 
locity listed in the NASA Extragalactic Database (NED). The 
latter was first transformed to the velocity relative to the lo- 
cal group (LG) centroid using the NED velocity calculator, 
and subsequently corrected for a LG infall towards the Virgo 
Cluster center {I = 284°, b = 74°) at 200 kms^^ (Tammann & 
Sandage a985< . 



2.2. Data acquisition 

The NIR images were observed with the ESO 3.6m NTT tele- 
scope. La Silla/Chile, during three consecutive nights from 
April 21st to 24th, 2000. The seeing ranged from 0'.'4 to I'.'O 
FWHM and the transparency variations were ^1% during 
the first and third night, and < 2% during the second night. 
Additional NIR images were taken at the 3.6m telescope of the 
German-Spanish Astronomical Center, Calar Alto, Spain, dur- 
ing several observing runs. The seeing conditions were gener- 
ally average (December 26th, 1999: FWHM 1'.'3 - 3'.'5, trans- 
parency fair; May 10th-15th, 2000: FWHM 0'.'8 - 1'.'4, trans- 
parency good to average; October 6th-10th, 2000: FWHM 
1'.'2 - 2'.'0, transparency good to fair). Both cameras used, the 
SOFI at the NTT Nasmyth focus and the OMEGA PRIME at 
the Calar Alto 3.6m telescope's prime focus, were equipped 
with 1024x 1024 pixel Rockwell HAWAII detectors. The pixel 
scales for SOFI, using the large field objective, and OMEGA 
PRIME were 0'.'292 and 0'.'396, yielding a FOV of 4'.94 and 
6'.76, respectively. The data was taken through the J and H 
broad band filters, as well as the modified K filters Kg at the 
NTT and K' at Calar Alto, both selected to extenuate the con- 
tribution of thermal background. 

To achieve an adequate sampling of the background vari- 
ations, the telescope was offset between exposures typically 
each 60 seconds. While compact galaxies were jittered within 
the FOV, larger objects required to alternately observe the ob- 
ject and the sky. To avoid detector saturation, images were ob- 
tained by summing up subexposures of few seconds each. The 
total on-object exposure times, after rejection of subexposures 
affected by unstable readout electronics or strong background 
gradients, are listed in Tabled 

2.3. Data reduction 

The exposures have been processed semi-interactively, em- 
ploying our NIR image reduction software based on ESO 
MIDAS ' . The procedures used for basic reduction, background 
subtraction, image alignment and coaddition follow the recipes 
which are detailed in, e.g., the SOFI User Manual^. 

Since surface photometry in the LSB regime depends sen- 
sitively on the quality of the background subtraction, care has 
been exercised to eliminate residuals of this correction on both, 
small and large spatial scales. For this purpose we have imple- 
mented a number of additional corrective steps into our soft- 
ware package. For each single science frame of an exposure se- 
quence we first computed an individual master background im- 
age, using subexposures that were taken closest in time, but at a 
sufficiently large distance from the target. Prior to the calcula- 
tion of the master background frame, the input background im- 
ages were cleaned from bright contaminating sources and nor- 
malized to the same mean intensity (see the SOFI User Manual 
for details). The master background frame computed this way 
was in turn scaled to the background level of the respective sci- 
ence frame, and subtracted. Each of the resulting subexposures 

' Munich Image data Analysis System, provided by the European 
Southern Observatory (ESO). 
" Issue 1 .2, available online at |http : / /www . eso . org| 
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Table 1. Sample galaxies 



Object 


RA(J2000) 




tj.NTT 


tff.NTT 


tXs.NTT 


Ms 


Ab 


D 


FWHM 


other 




DEC(J2000) 




tj.CA 


tH.CA 


tK'.CA 


(ref.) 




(ref.) 


(final)* 


names 








[s] 


[s] 


[s] 


[mag] 


[mag] 


[Mpc] 


["] 






(2) 




(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


Tol 3 


10''06'"33'= 


.6 


1440 


1560 


2100 


-18.0 


0.33 


13.8 


1.2 


NGC3125; ESO435-G041; 


(iE) 


-29°56'09' 


.0 


— 


— 


— 


(a) 




(a) 




Toll004-296 


Haro 14 


00''45"M6' 


.4 


— 


— 




-17.0 


0.09 


12.5 


1.6 


NGC 0244; UGCAOlO; 


(iE/nE) 


-15°35'49' 


.0 


1440 


1680 


1920 


(a) 




(a) 




VV728 


UM461 


ll''51'"33' 


.0 


1500 


1500 


2100 


-14.9 


0.08 


14.3 


0.7 


SCHG 1 148-020 


(il) 


-02°22'23' 


.0 


~ 


— 


— 


(i) 




(e) 






He 2-10 


08''36'"15' 


.0 


1140 


1500 


2040 


-18.7 


0.48 


8.7 


0.9 


ESO495-G021 


(iE/nE) 


-26° 24' 34' 


.0 


— 


— 


— 


(k) 




(f) 






Tol 1400-411 


14''03'"2r 


.0 


1140 


1500 


2100 


-16.5 


0.30 


4.8 


0.9 


NGC 5408; Tol 116; 


(iI,C) 


-41°22'44' 


.0 


— 


— 


— 


(1) 




(c) 




ESQ 325-G047 


Pox 4/Pox 4B 


ll'^51'"ll'= 


.6 


1500 


1500 


2100 


-18.8 


0.17 


46.7 


0.6 


CAM 1148-2020 


(il + comp.) 


-20° 36' 02' 


.0 


— 


— 


— 


(m) 




(b) 






Tol 65 


12''25"'46' 


.9 


2100 


2100 


3000 


-15.3 


0.32 


34.2 


0.6 


ESO380-G027; 


(iO/iI,C) 


-36°14'0l' 


.Os 


— 


— 


— 


(n) 




(b) 




Tol 1223-359 


Tol 1214-277 


12'M7'"17'= 


.1 


1200 


1500 


2100 


-16.9 


0.28 


102.6 


0.9 


SCHG 1214-277; 


(iO/iI,C) 


-28°02'33' 


.0 


— 


— 


— 


(0) 




(b) 




Tol 21 


Mkn 178 


11 ''33"" 29' 


.1 


240 


1380 


1380 


-13.9 


0.08 


4.2 


1.6 


UGC 06541; HOLMBERG 263 A; 


(iE) 


+49° 14' 17' 


.0 


— 


— 


— 


(P) 




(h) 




CGCG 242-046; SBS 1 130-1-495 


Mkn 1329 


12''37'"03' 


.0 


300 


360 


480 


-16.8 


0.10 


16.0 


0.7 


VCC 1699; IC 3589/91; 


(iI,C) 


+06°55'36' 


.0 


1200 


360 


960 


(i) 




(i) 




UGC 7790 


IC 4662 


17''47™06' 


.4 


1980 


1620 


660 


-14.9 


0.30 


2.0 


0.8 


ESQ 102-G014; 


(dl/iE) 


-64°38'25' 


.0 








(1) 




(g) 




He 2-269 


UM448 


11*" 42^12'= 


.4 


300 


360 


480 


-19.8 


0.11 


76.1 


0.7 


UGC 06665; SCHG 1 139-1-006; 


(il/BCG) 


+00°20'03' 


.0 


840 


900 


1200 


(1) 




(d) 




Ark312; ARP 161; Mkn 1304 



4|k: Resolution of the best image set available for the respective galaxy, after reduction and combination (a): Marlowe et al. 1 1999 1, (b): inferred 
from the heliocentric velocity t^hci listed in the NED, corrected for solar motion with respect to the center of the Virgo Cluster and adopting 
Ho=75kms"^Mpc"^ (cf. Sect.ITT}- (c): Karachentsev et al. <2002l, (d): Mirabel & Sanders 1988" (e): vlg from Smoker et al. ^20001. 
corrected for Local Group infall to the Virgo Cluster like b); (f): Tully <1988> . (g): Heydari-Malayeri 1990 (h): Schulte-Ladbeck et al. ( 2000>. 
(i): Yasuda et al. <1997> . (j): ms from Mendez & Esteban ( 2000^, (k): ms from Papaderos ( 1998 1, (1): ms from the RC3 (de Vaucouleurs et al. 
I1991> . (m): ms for Pox 4 from Mendez & Esteban U999.) ; no published ms for Pox 4B (n): ms from .P99. (o): ms from iFOl. (p): ms from 
Papaderos et al. <2002> 



was subsequently checked for a possible residual background 
gradient and, whenever necessary, rectified using a first order 
polynomial fit. The result image was obtained through combi- 
nation of all background-corrected and coregistered subexpo- 
sures. In this procedure, out of all pixels with identical sky co- 
ordinates, only those at sound detector coordinates have been 
included in the calculations. 

Small residual background variations in the final images 
needed to be interactively corrected by fitting twodimensional 
polynomials. Alternatively, regions still affected by small- 
scale residuals in the background subtraction were extracted 
as subimages, corrected in the latter manner, and inserted back 
into the original frame. Likewise, bleeding artifacts, caused by 
bright sources in the NTT/SOFI FOV, or blooming of bright, 
overexposed foreground stars in the very outskirts of IC 4662 
(northern edge) and Tol 1400-411 (east and west) were mod- 
eled and subtracted out. SBPs of the latter two galaxies were 
only analyzed well above intensity levels where slight residu- 
als from the replaced bright stars could possibly contribute. 



Final images taken in the same filter during different nights 
or at different telescopes were aligned to each other, trans- 
formed to equal pixel scales and resolutions, and coadded after 
weighting each one by its (S'/iV)^. The FWHM of the result 
images is listed in Tabled Images used for aperture or surface 
photometry were manually cleaned for fore- and background 
sources. 

2.4. Flux Calibration 

The SOFI data was calibrated by observing at different air- 
masses standard stars from Persson et al. (19981 six times each 
night. The excellent photometric stability throughout the NTT 
observing run has allowed for the derivation of an airmass- 
dependent calibration with a scatter of ^ 0.01 mag during 
nights 1 and 3 and ^ 0.02 mag during night 2. Zero points 
and airmass-dependent calibration coefficients both agree well 
with the average values supplied by the NTT/SOFI support 
team. Images of UM448 and IVIkn 1329, obtained by combin- 



Noeske et al.: NIR observations of Blue Compact Dwarf Galaxies 



5 



ing Calar Alto and NTT data, were calibrated using aperture 
photometry of bright stars in the SOFI FOV. Despite different 
K filters used at those telescopes, the integral fluxes are es- 
sentially preserved, as color terms in the filter transformations 
(Eqs.|3jl are not exceeding a few 0.01 mag. 

The BCDs Haro 14 and Mkn 178, for which no NTT ob- 
servations are available, were calibrated using the Two Micron 
All Sky Survey (2MASS) catalogue^ (Cutri et al. l2000l Jarrett 
et al. l2000l l. As 2MASS data (cf. Andreon l2002l may, due to 
their limited sensitivity, slightly underestimate the flux within 
the extended LSB component, we computed calibration terms 
using 2MASS field stars in the close vicinity of Haro 14 and 
Mkn 178. 

2.4.1. Transformation to other NIR photometric 
systems 

Unless stated otherwise, all magnitudes and colors given in this 
paper refer to the calibrations described in the previous sec- 
tion, and to the photometric systems defined by the instrumen- 
tal setup of each telescope (see Sect. |2}. However, whenever 
photometric quantities of the sample BCDs are compared either 
among each other, or with model predictions and data from the 
literature, they are first transformed to the 2MASS photometric 
system, using the relations described below. 

The calibration obtained at the NTT with the SOFI in- 
strument is based on standards by Persson et al. 1 1998|f, and 
can therefore be transformed to the Persson et al. < 19981 1 
Las Campanas Observatory (LCO) system using the color 
transformations given by the SOFI user manual (Issue 1.3, 
16/08/2000). A transformation from the latter system to the 
2MASS system is described by Carpenter ( 2001 1. Combination 
of both latter transformations yields the relations: 

i^s,2M =if.,s + 0.021(J -X)- 0.010 (1) 

(J - H)2M = 0.995( J - H)s + 0.015( J - K) + 0.002 

[H - Ks)2M = 1.029(i? - K,)s - 0.046( J - K) + 0.005, 

where the indices "2M" and "S" denote the 2MASS and SOFI 
systems, respectively. Cumulative uncertainties introduced by 
this transformation are <, 0.01 mag for Kg and ^ 0.015 mag 
for (J — H) and [H — Kg), given the typical range of values 
expected for NIR colors of BCDs. 

No transformations to other standard NIR systems are 
available for the photometric system defined by the Calar Alto 
3.6m telescope and the OMEGA PRIME camera. As stated in 
Sect. 12.41 the calibration of the Calar Alto frames is tied to the 
2MASS zero points. The remaining uncertainties introduced 
by the unknown color terms can be estimated from Carpenter 
(IMDtobe 

|(color)2M - (color)jyl < |0.08(color)o' | (2) 

per transformed color, i.e., the approximate upper limit of the 
color terms entering into the transformation between standard 
NIR photometric systems. 

^ [http : / / www . ipac . caltech . edu/ 2mass/| 



The K' and Kg magnitudes are related by the following 
equations (cf. Wainscoat Si. Cowie 119921 and the SOFI user 
manual): 

Kg^ K + 0.005( J - K) (3) 
K = K' - (0.22 ± Qm){H - K) 

2.4.2. Extinction correction: 

Magnitudes and colors given in this paper are corrected for 
Galactic extinction, adopting values derived from the B band 
extinction maps by Schlegel et al. ( 1998 1 (cf. Table[0 and the 
standard (Ry = 3.1) extinction law (Cardelli et al. 1989) im- 
plemented into the NED. No attempt was made to correct for 
internal extinction, since this is known to vary spatially even 
in the most metal-deficient BCDs (cf., e.g., Guseva et al. 1200 II 
Cannon et al. 120021 Hunt et al. l2003> and can be reliably con- 
strained in the SF regions only. 

3. Data analysis and SBP derivation 

3. 1. Derivation of surface brigtitness and color 
profiles 

Surface photometry aims at a standardized one-dimensional 
representation of a galaxy's two-dimensional flux pattern. One 
technique to compute surface brightness profiles (SBPs) re- 
quires the determination of the size A{pi) of the galaxy in □" 
for a series of surface brightness levels /i (mag/D"). By this def- 
inition, the equivalent radius R*=^J (A(/i)/7r) is a mono tonic 
function of the surface brightness and vice versa. In order to 
derive SBPs this way, one has to keep track of the morphol- 
ogy and angular extent of a BCD throughout its intensity span, 
i.e. in general to be able to interpolate an isophote down to the 
faintest measured level /i of an SBP. By this condition one can 
visually check for and screen-out fore- or background sources 
in the periphery of the galaxy, thus make sure that source con- 
fusion does not affect SBPs at faint levels. This task is more 
difficult to achieve when computing SBPs employing photon 
statistics inside circular or elliptical annuli, extending out to a 
user-defined maximal radius rmax- SBPs derived for an irreg- 
ular system by such techniques may vary from case to case, 
depending on, e.g., the adopted r,nax or "center" of the galaxy. 
It should be borne in mind that techniques of this kind, when 
applied to galaxies with a morphology significantly departing 
from the assumed circular symmetry, can strongly overestimate 
the exponential scale length and underestimate the central in- 
tensity of the LSB component (cf., e.g., Marlowe et al. ll997l l. 

In the present analysis, we compute SBPs using method iv), 
described in Papaderos et al. (2002 1. This is a hybrid technique, 
incorporating features of both aforementioned approaches (de- 
termination of the R* corresponding to a user-defined /i, as 
opposed to the determination of the mean surface brightness Jl 
inside a circular annulus with a user-defined radius i?*±Ai?*). 

For a set of n intensity intervals /„, with a mean intensity 
decreasing as n increases, masks Af„ are generated, each of 
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which extracts from a smoothed image of the galaxy the ar- 
eas with intensities within /„. The equivalent radius i?* corre- 
sponding to the mask area A„ is: 



Rl 



n— 1 -j^ 



(4) 



The surface brightness is calculated from the mean in- 

tensity of the original image within the area M„. Method iv) 
overcomes the problem of the artificial SBP flattening for low 
S/N levels (see discussion in e.g. Noeske 1999 Noeske et al. 
l2000l and Cair6s|2000i and proves reliable down to very faint 
surface brightness levels (/„ < sky noise). 

As a check for consistency, we also derived SBPs through 
ellipse fitting to isophotes or methods ii) and iii) in P96a apply- 
ing in all cases a minimum of image filtering to moderate pho- 
ton noise. The resulting SBPs have been compared with deep 
optical surface photometry to ensure that the underlying LSB 
component has been detected and modelled (see Sect. 13. 2> over 
a sufficient radius range. 

Color profiles were derived by subtracting the SBPs from 
each other, after the latter had been smoothed to equal FWHM. 
Surface brightness and color profiles, corrected for Galactic ex- 
tinction ('Sect. lZ4!2t . are shown in Figs.l2lthrough [T4l 

3.1 .1 . Notes on the errors 

The uncertainties of each data point include the fully propa- 
gated effects of (i) Poisson photon noise and (ii) residual back- 
ground variations on various spatial scales. The latter originate 
both from small residuals from the background reduction, and 
from undetected fore- and background sources. The treatment 
of (i) is explained in, e.g., P96a or Cairos et al. .2001a. The 
resulting uncertainties are in most cases small, owing to the 
large number of pixels typically involved in the calculation of 
each profile point at small S/N levels, i.e. in the outer LSB 
region of a galaxy. In these low S/N regions, the uncertain- 
ties from (ii) are typically dominant: background variations on 
spatial scales which are not small compared to the size of the 
considered mask segment or aperture do not average out. 

The amplitude of such variations was estimated as the stan- 
dard deviation of the mean intensities measured in quadratic 
apertures, placed on the background well away from the object, 
with sizes A of the order 1/10 of the object's diameter. The re- 
sulting uncertainties of each SBP point might be decreased by a 
factor of \J An/}?, accounting for the averaging of background 
variations over the area An (Eq.|4} of the respective mask seg- 
ment. This first-order correction was however discarded, since 
the background variations cannot be determined at the position 
of the galaxy's LSB component, and since the outer mask seg- 
ments cut through large areas of the image, possibly sampling 
individual stronger background anomalies. Error bars shown 
in the SBPs and color profiles (Figs. l2l through [141 therefore 
represent upper limits, i.e. typically overestimates, to the true 
errors, which are difficult to quantify. Nevertheless, the error 
constraints shown here give a better estimate of the true errors 
than the, typically too small, pure Poisson noise errors. 



The reliability limits of the SBPs were not estimated from 
these error constraints, but determined at a surface brightness 
where an isophote could still be robustly interpolated on a 
mildly smoothed image. 

3.2. Profile Decomposition 

The stellar emission of a BCD is due to the superposition of two 
distinct populations with respect to their A4/C ratio and spatial 
extent: (i) the underlying LSB host galaxy which, owing to its 
high Ai /C, contains the bulk of the system's stellar mass, and 
(ii) the younger stellar population, attributable to the ongoing 
and recent SF activity. In the main class of iE/nE BCDs, the 
younger stellar component dominates the optical emission in- 
side ~2 exponential scale lengths of the underlying LSB host 
and contributes, together with ionized gas emission, ~l/2 of 
the B light within the 25 B mag/D" isophote (P96bl Noeske 
1999 Salzer & Norton fT999V Evidently, a meaningful study of 
the structural properties of a BCD requires the decomposition 
of SBPs into these two main photometric components. In the 
following, we employ a simple SBP decomposition scheme, in 
which we fit only the LSB emission. Subtraction of the best fit 
LSB model from the SBP allows us to deduce the luminosity 
fraction and spatial extent of the superimposed SF component. 
In order to make sure that extended starburst emission does 
not affect the decomposition results, we fitted the LSB com- 
ponent beyond a transition radius (cf. IP96al IC03I . where 
optical and optical-NIR color gradients vanish and isophotes 
become more regular In addition, whenever available, we used 
Ha maps to trace the size of the SF component, and get from it 
an additional constraint to R^. The LSB emission has been fit 
out to the radius where SBPs became increasingly uncertain as 
a result of background noise (typically for ^O.Subgr for method 
iv). The fitted SBP range in J is indicated with the dashed-gray 
line at the bottom of each profile (Figs.l2HT4l. 
An exponential fitting law 



I{R*) = /e,o exp 



R^ 

a 



(5) 



with a central intensity /e.o and an exponential scale length a 
has been found to approximate well the intensity distribution of 
our sample galaxies in their LSB component. An exponential 
law is also known to fit well the outer parts of dwarf uregulars 
(dis, Patterson & Thuan [T996l van Zee|2p00l and dwarf ellipti- 
cals (dEs, Vigroux et al. ll988l Binggeli 8c Cameron ll991l l over 
a few scale lengths. 

However, for several of our sample BCDs, an extrapola- 
tion of the exponential LSB slope to smaller radii yields for 
i?*^l. . . 3q; a higher intensity than the observed value. Thus, 
a meaningful decomposition of such SBPs cannot be achieved 
on the usual assumption that the exponential law is valid in 
the LSB component all the way to R*=Q". Instead, one has to 
adopt an alternative fitting formula, approaching the exponen- 
tial law for large radii and flattening in its inner part. Such a 
distribution should be compatible to the "type II" profiles of 
disc galaxies (Freeman 1 19701 see also MacArthur et al. l2002> . 
or the "type V" profiles described for spheroidal early type DGs 
by Binggeli & Cameron (tl991j . 
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Name Band fiE,o a mLgB Piso JTiPiao Eiso niEj.o msBP mtot rcffsrso ?7sbp 

(b,q)'' mag/D" kpc mag kpc mag kpc mag mag mag kpc 

0) (2) (3) (4) (5) (6) (7) (8) (9) (10) (11) (12) (13) 

TolS" J 18.24±0.05 0.52±0.01 11.78 1.20 12.42 2.29 11.86 11.31 11.29 0.61,1.23 1.45 

■k H 17.74±0.06 0.53±0.01 11.26 1.18 11.76 2.07 11.37 10.73 10.69 0.61,1.22 1.43 

Ks 17.45±0.09 0.51±0.01 11.06 1.18 11.48 2.12 11.15 10.50 10.46 0.57,1.15 1.49 

Haro W J 17.46±0.12 0.37±0.01 12.70 L39 12.45 TSS 12.82 TlSI 11.78 0.61,1.21 TTs" 

3.6,0.94 H 16.97±0.26 0.39±0.02 12.11 1.37 11.69 1.76 12.31 11.12 11.07 0.61,1.22 1.25 

0.49 K' 16.88±0.94 0.39±0.08 11.99 1.32 11.54 1.83 12.17 10.98 10.97 0.60,1.21 1.45 

UM46r J 19.37±0.18 0.21±0.01 15.56 042 16.19 069 15.88 15.09 15.04 0.34,0.57 r96" 

2.3,0.85 H 18.82±025 O21±0.02 14.98 0.35 16.12 0.60 15.46 14.66 14.57 0.37,061 1.82 

0.49 Ks 18.47±079 O19±O04 14.79 036 15.75 063 15.14 14.46 14.37 0.32,054 2.04 

Henize2-10'^ J 18.46±0.03 O67±O01 10.46 087 1057 Z80 1055 9m 9?70 0.41,1.26 ISO" 

★ H 17.75±O04 0.64±0.01 9.85 0.78 9.94 2.51 9.96 9.17 9.13 0.4O1.22 2.44 

Ks 17.51±O03 O62±0.01 9.68 0.79 9.62 2.56 9.78 8.95 8.89 034,1.04 2.44 

Tol 1400^1 P J 18.29±011 O33±O01 11.18 069 12.67 Tm 11.32 1099 1091 0.60,0.98 iM 

3.0,0.82 17.30±O14 O27±0.01 1057 0.31 13.65 1.18 1071 1057 1055 0.57,093 1.33 

063 Ks 17.39±0.09 O30±O01 10.43 O50 12.49 1.28 1058 10.32 1031 059,097 1.75 

Pox 4" J 19.26±015 0.86±0.03 15.15 Z27 14.80 IST 15.54 14.24 14.20 1.24,2.17 OF 

2.7,0.90 19.10±O31 O88±O07 14.95 2.18 14.35 1.79 16.10 13.85 13.81 1.27,2.15 1.45 

0.42 K/ — — _____ 13.54 13.52 1.18,2.09 1.96 

Pox4B*= J 19.47±0.40 0.29±0.04 17.37 OSl 19.29 093 17.67 VfTl 17.12 0.54,0.87 0.75 

2.0,0.80 U"^ — — _____ 16.83 16.72 0.51,080 0.77 

0.49 K/ — — — — — — — 16.75 16.52 0.47,069 0.75 

Tol65*= J 19.58±018 0.28±0.02 17.25 0.43 17.47 0.82 17.79 16.57 16.55 0.37,080 1.59 

2.7,0.90 H'^ — — _____ 16.23 16.30 0.42,084 — 

0.46 Jf/ — — — — — — — 16.03 15.99 0.24,061 — 

Tol 1214-277'=''' P 19.16±024 O53±0.03 18.05 075 18.56 1^69 18.55 17.46 1741 0.86,1.83 r89" 

3.3,0.92 H'* — — — — — — — 17.22 17.26 1.15,1.90 — 

0.48 K/ — — _____ 16.94 17.07 0.81,1.67 — 

Mkn 178-' J 2062±011 0.28±0.01 12.94 039 14.36 062 ViM 12.78 12.70 0.34,0.57 r27" 

★ H"* 19.87±018 O26±0.02 12.35 0.26 14.85 0.51 12.92 12.29 12.25 0.34,058 1.12 
K"* 19.92±0.43 O27±O04 12.31 036 14.00 0.52 12.91 12.13 12.01 036,058 097 

Mkn 1329'= J 19.25±O04 O67±O01 \2M 099 1461 232 iJoi 12.64 12.60 1.12,1.97 ToT 

1.6,0.70 H'* 18.81±O01 O71±0.01 12.25 1.01 13.85 2.09 12.57 12.00 11.94 1.17,2.08 1.11 

068 Ks 18.24±0.06 O60±O01 12.03 079 14.56 2.09 12.23 11.95 11.93 1.09,1.87 072 

IC 4662'= J 16.79±024 0.15±0.01 1O60 055 1032 083 1073 9^68 9^64 0.29,053 lAl 

4.8,0.97 H 16.61±073 O15±0.02 1036 0.58 9.64 069 1066 9.15 9.05 0.29,052 1.37 

0.48 — — _____ 9.04 8.93 0.27,0.48 1.75 

TJM44F J 20.55±0.07 3.83±0.10 13.48 426 12.50 8^66 13.92 12.17 12.15 1.28,3.63 233" 

★ H 19.95±017 4.02±0.31 12.77 4.08 11.88 7.60 13.39 11.52 11.50 1.32,3.94 2.40 
Ks 19.55±0.15 3.61±019 12.61 4.01 11.53 8.16 13.05 11.22 11.30 1.21,3.24 2.44 



a: All values are c orrected for Galactic extinction, adopting the An from TablellKcf. Sect. l2.4.2t . 

b,c: See Sect. l3.2i for details. Objects whose LSB component has been modelled by a pure exponential fit (Eq.|5} are marked with an 
asterisk. 

d: Decomposition uncertainty higher (faint object, or data affected by nearby bright stars). 

e: Calibration obtained at the ESO NTT, cf. Sect.l2n 

/: Calibrated using field stars from the 2MASS catalogue, cf. Sect. l2.4l 

g: Fit radius interval and (b,q) adopted from lFOll 



Previous optical studies have allowed for the detection and 
modelling of such an inwards flattening exponential LSB pro- 
file in a few BCDs only (e.g. I Zw 1 15. IP96al Tol 65, Papaderos 
et al. 1999, hereafter 'P99', Tol 1214-277, Fricke et al. 2001, 
hereafter FOl SBS 0940H-544, Guseva et al. 12001 ). However, 
SBPs of this kind do not appear to be rare among intrinsically 
faint dEs and dis. They have been observed in several dEs with 
and without a central nucleus (cf. Binggeli & Cameron . 1 99 !■ 



Cellone et al.T^ Young & Currie'1994"1995'), and Vennik et 
al. (.20 001 deduce a fraction of >10% for the SBPs of late-type 
dwarf galaxies falling into this category. In fact, a casual in- 
spection of the sample of Patterson & Thuan ( 19961, Makarova 
et al. (19981, Vennik et al. (2000 1 and van Zee (2000) reveals 
several examples of dIs showing an exponential outer LSB 
slope and a pronounced flattening for intermediate to small 
radii (e.g., UGC 2034, UGC 2053, UGC 5423, UGC 9128, 
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UGC 10669 in the sample of Patte rson & Thuan lT996l KKH35 
and KKSG 19 in Makarova et al. [T998l or UGCA 9, UGCA 
15, UGC 2345,UGC 9240, UGC 10445 in the sample by van 
Zee 12000 ). Other examples are the dis Holmberg I (Ott et al. 
l200n , Holmberg II (Noeske |T999j, H 1032-2722 (Due et al. 
11999 1 and Kar 50 (Davidge 2002 1. Following the nomencla- 
ture of Binggeli & Cameron Cl991j we shall henceforth refer 
to this type of SBPs as to type V. 

The physical origin and exact form of typeV profiles in 
puffed-up stellar systems has not yet been studied in detail, nei- 
ther observationally nor theoretically. Therefore it is difficult 
to say which fitting formula approximates best their intensity 
distribution (see discussion in Sect. I5.3l l. An empirical fitting 
function that yields by deprojection a finite luminosity density 
for R*=Q" has been proposed in P96a, as 

I{R*) = /E,oCxp (^-^^ {1 - gexp(-P3(i?*))} (6) 
with 

The modified exponential distribution Eq. (|6j, in the following 
referred to as med, flattens with respect to a pure exponential 
law inside of a cutoff radius ba, and attains at R*=Q" an inten- 
sity given by the relative depression parameter q — A///e,o < 
1. An advantage of the med is that its exponential part (left part 
of Eq.|6}, depending on /e,o and a only, can be separately con- 
strained by fitting Eq.|5lto the outer exponential part of a type V 
profile. This has been done by means of unweighted fits, after 
SBPs were splined to equidistant radius steps, in order to ensure 
that the fit solution is not biased through a clustering of points 
to a specific radius range. Once the scale length a is fixed, the 
multiplicative right-hand side part of Eq. (|6|l is a function of 
b and q, only. Because the {b,q) parameter space has not been 
explored so far, e.g., by fitting Eq. (|6j to a sample of dEs with 
type V SBPs, it is difficult to judge which values are appropri- 
ate for dwarf galaxies falling into this category. P96a inferred 
a tentative ratio of b/q ~3, they noted, however, that this could 
be luminosity-dependent. 

In order to deduce plausible constraints to {b,q), we first 
subtracted from the J image of each sample BCD most of the 
irregular starburst emission. The latter has been approximated 
by an unsharp-mask version of the original image (Sect. I3.4> . 
adjusted such that no regular emission from the LSB compo- 
nent was removed. Ha exposures and color maps were used 
to further ascertain that the subtracted emission was not part 
of the underlying host galaxy. J band exposures processed this 
way were then used to compute SBPs for the LSB component. 
These profiles, denoted Jlsb, allow to better trace the LSB 
component down to smaller R* than the ones derived prior to 
partial subtraction of the starburst light. By fitting Eq. (|6} to 
the Jlsb SBPs we derived b,q (first column in Table|2}. In most 
cases, fit uncertainties are 0. 1 and 0.05 for h and q, respectively. 
No reliable Jlsb profiles could be computed for Tol 1214-277. 
For this system we fixed {b,q) to values inferred bv lFOll from 
optical VLT data. 



Alternatively, typeV Jlsb profiles were fit with a Sersic 
model (Sersic 1968 1 of the form 

/(i?*) =/s,oexpU— j . (8) 

The exponent^ 77 in Eq. (|8} is a handy indicator of systematic 
deviations of the observed SBP from an exponential distribu- 
tion. An 77 < 1 corresponds to a convex profile, which might be 
compatible to a med distribution with q> Q (see discussion in 
Sect. 15. 3> . whereas an 77 > 1 indicates a concave profile which, 
in the special case of 77 = 4, translates into the de Vaucouleurs 
law (see, e.g., Caon et al. 1993 1. 

Note that Eq. (|8} approximates well the projected light of a 
variety of stellar populations with both a nearly constant M. / C 
(ellipticals or bulges; Caon et al. 1993 Andredakis et al. 1 19951 
Graham et al. 119961 early- type dwarfs; Cellone et al. 119941 
Young & Currie 1994) or strongly varying / C (for instance, 
nE BCDs or the plateau component in the SBP of iE BCDs; 
IP96al l. Therefore a Sersic exponent 77 > 1 gives no strong 
indication for a stellar system being similar to an intrinsically 
luminous early-type galaxy (see also discussion in C03 ). 

The Sersic exponents of the LSB components (77LSB) were 
obtained from non-weighted fits to the Jlsb SBPs (see above) 
of type V profiles. This limits 7/lsb to values <1. We avoided 
to fit Eq. (jSJl to the outer part of SBPs alone, i.e. for R*>R^, 
as solutions obtained this way are very uncertain and depend 
strongly on the accuracy of the sky subtraction (see Sect. 15.31 
and detailed discussion in lC03t . Notwithstanding the fact that 
Eqs. (|6j and (|8j give comparably good fits in terms of x^, we 
decided not to include the full {Is,o,P,il) Sersic solutions for the 
Jlsb SBPs in Tabled Instead, in its Col. 1, we quote only the 
Sersic exponent 77LSB, the uncertainties of which are estimated 
to be of the order of 20%. 

The photometric quantities of the sample BCDs are sum- 
marized in Table|2l BCDs without signatures of a type V profile 
in their underlying LSB component are marked with an aster- 
isk in column 1. For the remaining systems we list the (b,q) 
and 7/LSB parameters, as obtained respectively by fitting Eqs. 
(|6j and ^ to Jlsb SBPs. Columns 3 and 4 list, respectively, 
the extrapolated central surface brightness /iE,o (mag/D") and 
exponential scale length (kpc), obtained by fitting Eq. (|5} to 
the outer exponential LSB part of each SBP. Column 5 lists the 
total apparent magnitude of the LSB component, computed by 
extrapolating the fitted model (i.e. Eq.|5lor Eq.|6j to R* — 00. 
Columns 6 through 9 list the radii and magnitudes of the star- 
forming (P) and underlying stellar LSB component (E), as ob- 
tained by profile decomposition. Following P96a, we measure 
the respective radial extent (Piso, ^'iso) and encircled magni- 
tude ('mp^^^,mEis„) of each component at an isophotal level 
iso, taken to be 23 mag/D" for J and 22 mag/D" for H and 
K. The isophotal radii determined for the sample BCDs at 23 
J mag/D" turn out to be comparable to those obtained from 
optical SBPs at 25 B mag/D" (P25 and E25 in P96a). Column 

■* For consistency with, e.g., Caon et al. il993l and iC03i we define 
here the Sersic shape parameter as rj, i.e. the reciprocal exponent in 
Eq. Jsj. Note that in other studies (Cellone et al. ll994l Young & Currie 
IT994..1995..P96aj the Sersic shape parameter is referred to as I/77. 
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10 lists the magnitude from an SBP integration out to the last 
data point, and total magnitudes from aperture measurements 
(cf. Sect. 13. 3t are listed in column 11. The radii r^ff and rgo, 
enclosing 50% and 80% of the SBP's flux are included in col- 
umn 12. Finally, a formal Sersic exponent for the whole SBP 
(VSBp), for later comparison with literature data, is listed in 
column 13 of Table|2] 

Errors of the fitting law parameters determined from un- 
weighted fits are an underestimation of the true errors. LSB 
slope differences between J, H and K (Table |3 should there- 
fore not be considered significant, but rather to reflect the true 
exponential slope uncertainties, which are typically up to 10% 
in J, and may become somewhat larger in H and K, depending 
on the data quality. 

3.3. Aperture photometry and total magnitudes 

Since a SBP can only be accurately derived in regions with a 
sufficiently high S/N (i.e., in general, down to the minimum 
intensity level in which the morphology of a BCD can still be 
visually checked and the problem of source confusion can be 
handled; cf. Sect. 13. 1> . profile integration out to the last mea- 
sured data point may, in some cases, underestimate the object's 
total flux. For extended LSB sources, or when NIR SBPs do 
not go sufficiently deep (e.g., the SBP of Tol 1214-277), the 
fractional flux missed can easily exceed 10%. Similar problems 
may affect growth curve flux determinations, since these basi- 
cally rely on a crude SBP derivation and can sensitively depend 
on the quality of background subtraction. 

We therefore measured total magnitudes within polygonal 
apertures which extend typically out to 1.5 Holmberg radii 
(Col. 11 of Table |2li, after removal of fore- and background 
sources from the area of interest (cf. Sect. 12. 3> . Errors take into 
account the Poisson noise and small-scale variations of the lo- 
cal background. 

Magnitudes and colors of selected features, such as stellar 
clusters or H II regions, have been corrected for the flux contri- 
bution of the underlying LSB component, by interpolating the 
mean surface brightness of adjacent regions. Values computed 
this way are marked with the superscript if. As pointed out in 
Cairos et al. ('2002a) and Papaderos et al. 2002 corrections for 
the LSB emission are generally not negligible in optical wave- 
lengths. That this statement is also true in the NIR domain is 
illustrated on the example of the iE BCD UM 461 (Fig.^; cor- 
rection for the flux contribution of the LSB background shifts 
even the brightest SF region (a, mj|=17.26 mag, open circle) 
of this system by H-0.25 mag and -0.1 mag in the J — H and 
H — Ks color diagram. A census and photometric study of 
compact stellar clusters in the extended BCD sample included 
in this project after correction for the LSB- and ionized gas 
emission will be presented in a subsequent paper 

3.4. Unsharp masking 

The considerable intensity range of a BCD, from its faint LSB 
outskirts to the brightest nuclear starburst region, renders the 
detection of fine coherent morphological features in the central 
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Fig. 1. NIR two-color diagram for the brightest compact re- 
gions in the il BCD UM461. Filled circles show the colors 
of the regions a through g indicated in Fig.|5j after correction 
for the flux contribution of the underlying LSB host galaxy. 
Open circles, connected with dotted lines, indicate the color of 
the respective region prior to that correction. The star marks 
the color of the stellar host galaxy ("LSB"). The color range 
covered by red supergiants in the SMC (Elias et al. I1985t is 
shown by the box labelled RSG to the lower right. The tem- 
poral evolution of the J — H vs. H — Kg color for an instan- 
taneously formed stellar population with Zq/IQ for an age be- 
tween <=0 and t=50 Myr (solid line) has been calculated with 
the PEGASE code (Fioc & Rocca-Volmerange [T997t . The col- 
ors for t =0 correspond to purely gaseous emission. Arrows de- 
pict the effect of a 10% contribution of warm ('^200 K) dust to 
the K band (cf. Campbell & Terlevich 19841, as well as the ex- 
tinction vector All colors are transformed to the 2MASS pho- 
tometric system. 



portion of the galaxy difficult. We use therefore a modified un- 
sharp masking technique (cf., e.g., Papaderos 1 1998t . referred 
to in the following as hierarchical binning (hh) - transforma- 
tion. This contrast-enhancing procedure is stable against noise 
at low intensities and allows for the flux determination of faint 
sources within the bright background of a BCD, with an angu- 
lar size smaller than a user-defined value. Morphological fea- 
tures of interest revealed using this procedure are displayed in 
the grayscale/isophote insets of the sample galaxies, and de- 
scribed in Sect.|4] 

3.5. Colors of the underlying LSB host galaxy 

Colors of the LSB host galaxy were derived as the error- 
weighted mean of the color profiles for R*>R^, after rejec- 
tion of deviant points, being probably affected by uncertain- 
ties in the sky determination, and local residuals in the subtrac- 
tion of background sources. For Tol 1400-41 1, Pox 4, Mkn 178 
and IC 4662, uncertainties in the LSB colors are larger, due 
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Table 3. Colors of the host galaxy' 



Object 


J-H 








[mag] 


[mag] 


[mag] 


Tol3 


0.55±0.07 


0.10±0.06 




Toies^^ 


0.48±0.17'' 




1.3±0.15'' 


Tol 1214-277'= 


0.39±0.16'' 




1.0±0.13'' 


Tol 1400-411'= 


0.43±0.08 


0.26±0.09 


1.45±0.16 


Pox 4'= 


0.42±0.13 


0.12±0.13 




Pox4B 


0.43±0.15'' 


O.liO.15'' 




UM448 


0.69±0.10'' 


O.MzbO.U'' 


2.09±0.15 


UM461 


0.50±0.06 


0.15±0.06 


1.92±0.16 


He 2- 10 


0.59±0.06 


0.11 ±0.07 




1C4662'= 


0.57±0.12'' 


0.03±0.12'' 




Mknl78 


0.53±0.15 


0.25±0.16'' 


2.00±0.15 


Mkn 1329 


0.65±0.09 


0.03±0.10 




Haro 14 


0.69±0.12'' 


0.20±0.12'' 


2.02±0.15 


Corrected for galactic extinction (see Sect.l2.4.'2t: NIR colors are 



transformed to the 2MASS system (cf. Sect. 12.4. H . The B band 
denotes the Johnson B. Errors include cumulative uncertainties in the 
determination of color profiles, the calibration and the transformation 
to the 2MASS system. 

possible local instabilities in one SB? at low S/N levels. 
'=: possible contamination by gas emission over a large portion of the 
LSB host galaxy 



to extended starburst emission or crowding with nearby bright 
stars. As for the very metal-deficient systems Tol 65 and Tol 
1214-277, the faintness of their LSB component in Kg has not 
permitted, despite generous exposure times, to pin down their 
H — Kg colors. 

Whenever calibrated optical data were available, optical- 
NIR colors were derived. Because the quality of the SBPs was 
typically better in J than in H, and the J — H color shows 
little evolution with time (few 0.1 mag) for old (^ IGyr) stel- 
lar populations, we derived B — J colors instead of the more 
commonly used B ~ H colors. 

The mean colors of the host galaxies are shown at the right 
edge of each color profile (Figs.l2l-I14>. in the photometric sys- 
tem in which the respective galaxy was observed and calibrated 
(see Sect. l2.4.ll . Table|3llists the NIR colors, transformed to the 
2MASS photometric system to facilitate comparisons and the 
B — J color, where available. Errors give cumulative uncertain- 
ties in the calibration, transformation to the 2MASS system, 
and the scatter and systematic uncertainty of each color profile. 

4. Results and discussion of individual objects 

4.1. Tol3 (NGC3125) 

This luminous {Mb 18.0; Marlowe et al. [T999t and 

relatively metal -rich ( Zq/6. . . Zq/3; K obulnic ky et al. [19991 
Schaerer et al. 119991 Marlowe et al. "1999' and references 
therein) BCD is known to be a member of the NGC 3175 
galaxy group (Garcia [1993}. NIR and optical images reveal 
two compact {<, 0.3 kpc) high-surface brightness (HSB) re- 
gions: the brighter northwestern knot A, roughly coinciding 




£0 7.C 



Fig. 3. top: J — H color map of the central region of Tol 3, 
corrected for Galactic extinction. The overlaid J contours are 
computed from the contrast-enhanced blow-up included in Fig. 
[2|(see text for details), bottom: EW(Ha) map of Tol 3, overlaid 
with R contours. 



with the geometrical center of the smooth LSB host galaxy, and 
the fainter knot B, located ~1 kpc southeast of A (see Fig. [2^). 
Either knot is the locus of ongoing star formation, as witnessed 
by the detection of Wolf-Rayet features (Kunth & Sargent 
[Ml] Vacca & Conti [19921 Schaerer et aL [T999l and red su- 
pergiants from CO absorption studies (Campbell & Terlevich 
.1984j . The intense SF activity in Tol 3 is also reflected on rela- 
tively blue optical colors (B - V=0.24, V - 7=0.28; Marlowe 
et al. I1997> in its nuclear region, as well as on copious Ha 
emission (1.1 x 10"'^ erg s^^; Marlowe et al. 1997). Deep Ha 
imaging by Marlowe et al. ( I1995> revealed a bipolar outflow 
roughly perpendicular to the major axis of the BCD, extending 
out to ~2.8 kpc from its nuclear region. 

Unsharp-masked NIR images reveal a complex morphol- 
ogy in the HSB regime of the BCD (inset in Fig.|2t). Regions 
A and B are immersed in an extended "S-shaped" pattern, ^0.9 
kpc in length, ending at its NW tip with a curved feature (N 
arc). There is some evidence for propagation of SF activi- 
ties, as both archival optical NTT and NIR data reveal signa- 
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CrOSOC r' [arcsec] 

Fig. 2. a): Contours overlaid with a J image of Tol 3 (Z)=13.8 Mpc). North is up, east to the left. Contours, corrected for Galactic 
extinction, go from 19 to 23.5 J mag/D" in increments of 0.5 mag. The 23 J mag/D" isophote is illustrated by the thick contour. 
The brighter star-forming region A and the fainter knot B (following the denomination by Schaerer et al. 1 19991 1 are marked with 
crosses. The inset shows a contrast-enhanced (see Sect. 13 .4t and magnified version of the central region of the BCD (indicated 
by brackets in the contour map image). Compact sources arranged along the southwestern arc -like chain are marked with small 
crosses, b): Surface brightness profiles (SBPs) of Tol 3 in the J, H and Kg, corrected for galactic extinction. For a better 
visualization, the H and Kg SBPs are shifted by -1 and -2 mag, respectively. The thick solid line illustrates an exponential fit 
to the stellar LSB component in J (cf. Sect. 13. 2> . computed in the radius range indicated by the light gray, long-dashed bar at 
the bottom of the figure. The emission in excess to the fit (small open circles) is attributable to the starburst component, which 
dominates the light in the inner part of Tol 3. The isophotal radii P23 and £'23 of the star-forming and LSB component at the 
surface brightness level of 23 J mag/D" (horizontal dotted line) are indicated. The bar at the upper left of the figure corresponds 
to a galactocentric distance of 1 kpc. (c): Color profiles, computed by subtraction of the SBPs shown in the upper right panel. 
The thick lines at the rightmost part of the diagram indicate the mean J — H and H — Kg colors of the LSB component (see 
Sect.l33l. 



tures of a younger age and stronger ionized gas emission to- 
wards region B (cf. Fig.|3j bottom). For region A we determine 
within a rectangular 4"x4" aperture colors of J — H=0.6l 
( J - Hi=0.56) and H - K=OAl (H - K^=0.5^). For knot 
B we infer a J - H=0.39 (J - i/t=0.14) and H ~ K=OAS 
(H — K^=Q.ll). Such colors suggest a younger stellar age to- 
gether with an appreciable ionized gas contribution towards the 
SE part of the SF component. The latter is verified from the Ha 
equivalent width (EW) map in Fig. |3l (bottom) (see also Gil de 
Paz et al. 2002) which shows that SE of knot B, and all over an 
extended rim perpendicular to the major axis of the BCD, the 
EW(Ha) rises to >200 A. 

Interestingly, unsharp-masked NIR and optical images re- 
veal on larger scales a chain of faint knots (depicted with 
crosses in the inset of Fig.|2t), aiTanged over ^2.5 kpc SW of 
regions A and B. Their typical magnitudes of '^20.6. . . 19 
mag translate into absolute magnitudes of -10. . . -12 mag. The 
nature and formation history of this extended feature is in- 



triguing. One interpretation is that it delineates the approach- 
ing side of an oblate star-forming shell triggered by the burst, 
or that it may be associated with an inclined circumnuclear 
disk of ^1 kpc in radius. The first hypothesis is consistent 
with the findings by Alton et al. 11994^, who suggested from 
imaging-polarimetry the presence of a large-scale bipolar re- 
flection dust-nebula illuminated by the central starburst region. 
This could provide effective UV-shielding, thereby allowing for 
secondary SF activity. The systematically redder colors in the 
NE half of Tol 3 (Fig.|3j top) are also in line with the same hy- 
pothesis, if they originate from the more strongly absorbed far 
side of the galaxy. 

SBPs in the NIR (Fig.|2b) show in the radius range 22" < 
R* <48" an exponential intensity fall-off with a scale length 
a — 0.52 kpc. This value is close to the B scale length of 
a =0.48 kpc, inferred for the LSB component by Marlowe et 
al. f 1997) within ^30". By contrast, our surface photometry 
does not appear- to be compatible to that of Kunth et al. (.1988.1 . 
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These authors show optical SBPs out to a radius i?*=80", by 
a factor of 1.6 larger than the study here or in Doublier et al. 
( 1999 1 and up to 2.7 times larger than in Marlowe et al. ( 1997 1. 
NIR SBPs show no evidence for a dominant W-/^ profile in Tol 
3 (cf. Kunth et al. 1988 Doublier et al. 1999 1. A de Vaucouleurs 
profile is neither supported by the Sersic index 7y = 1.4 we 
derive for the entire J profile. 

The NIR color profiles (Fig. |2j;) show minor gradients 
(<0.15 mag kpc^^) and level off to J — _ff=0.54 mag and 
H — K=Q.\2 mag for R*> 20". These results are not compati- 
ble to Doublier et al. J2001I I. who report for j?*>10" a roughly 
linear J — H color increase from ^0.5 mag to ~2.0 mag. A 
good agreement is found with Vanzi et al. ("2002"), who derive 
for the LSB component of Tol 3 colors of J — H=Q.6 mag and 
H - if =0.25 mag. 

4.2. Ham 14 (NGC 244) 

This relatively metal-rich nE BCD {Z^Zq/3, Hunter & 
Hoffman 1999 1 shows SF activity ontop a smooth, nearly circu- 
lar stellar LSB host galaxy (Fig.|^). The intensity distribution 
of the latter (Fig.|4j3) is approximated best by a modified ex- 
ponential distribution (Eq. |6|l, flattening for R*^2Q" (see the 
detailed discussion in Sect. 15. 2> . 

The central part of the BCD contains a massive complex of 
SF regions, displaced by ~0.5 kpc SW of the geometrical cen- 
ter of the outer LSB isophotes. This SF component measures 
Ril.4 kpc in diameter in the J band (this paper) and ~2 kpc in 
the Ha line (Marlowe et al. ll997> . Its composition out of sev- 
eral individual regions, reported by Doublier et al. (1999 1 from 
optical images, is also evident from the contrast-enhanced NIR 
images (Fig.|^, upper-left inset), which reveal a wealth of in- 
dividual regions, spanning a range of ^3 mag. The upper-right 
inset of Fig. |^ reveals a chain of faint knots with a length 
of w 1.5 kpc (labeled "E loop"), extending eastwards from 
the main SF complex. There is a hint for a similar feature in 
the southwestern direction. Published data (e.g., the Ha image 
by Marlowe et al. Il997> do not allow to assess whether these 
faint features may trace induced star formation along super- 
giant shells, as might be hypothesized from their morphology. 

For the two brightest regions, denoted a and b, we obtain 
respective absolute magnitudes of Mj-'-=-14.7 mag and -13.3 
mag, and effective radii <80 pc. If the mean E{B — V) =0.35 
for Haro 14 (Hunter & Hoffman 1999 ) applies to knots a and b, 
then their de-reddened colors would be J — and H — 

Ks^ ^0-2. Such colors are reached at the earliest when the NIR 
emission becomes dominated by red supergiants (t'^lO. . .30 
Myr), and suggest no strong nebular line contamination. This 
points against substantial ongoing star formation in regions a 
and b. 

The large spatial extent of SF sources in the inner portion 
of the BCD is further evidenced by the profile decomposition 
(Fig. 13)), yielding a plateau radius P^«1.4 kpc in the J band. 
The optical and NIR colors of the LSB host galaxy indicate 
an age of several Gyr, in agreement with previous estimates by 
Marlowe et al. (il999<. 



4.3. UM461 

Telles & Terlevich ( fT995t suggested that the il BCD UM 461 
forms together witii UM 463, UM 465 and UM 462 a loose 
group of dwarf galaxies. 

SF activity is confined to the whole northeastern part of the 
galaxy, i.e. within the 21.5 J mag/D" isophote, or on a spa- 
tial scale of ^--^0.9x0.7 kpc. The two brightest SF regions a 
(mj*=17.3 mag; Fig.|5li) and C (mj*=19.5 mag) in UM461 are 
separated by 20 km s^^ in the velocity space. This difference 
is of the order of the intrinsic H I velocity dispersion within UM 
461 (-30 km s"^ van Zee et al.'T998't. The available data al- 
low us to resolve a manifold of morphological features within 
the SF component, most notably a chain of compact sources 
northwest of region C (labelled e-g) with a typical Mj^^ - 
10.4 mag. Their J — colors of 0.5-0.6 mag are consistent 
with the interpretation by Mendez & Esteban ( 2000 1 that they 
are extinguished SF regions, formed not earlier than 100 Myr 
ago. From profile decomposition we infer the summed up lu- 
minosity fraction of compact and diffuse SF sources to —37% 
of the J band light of UM 46 1 . 

UM 46 1 shows in its faint outskirts a slight asymmetry to- 
wards the SW direction. The intensity profile of the LSB host 
can be approximated by a med model with a scale length 
a=0.21 kpc and a depression parameter q «0.85 (Fig. (Sji). 
The J band scale length derived here is in excellent agree- 
ment with the value derived in the optical by Telles et al. 
fl997i. Color profiles reflect the ongoing SF for small radii 
(R*^2"), and approach mean values of J — H =0.49 mag and 
H — Ks =0.17 mag in the LSB periphery. Such colors, together 
with the B — J K,\.9 mag determined from optical data, point 
consistently to a relatively evolved stellar LSB background. 
The NIR colors derived here do not appear to be compatible 
with those by Doublier et al. (2001 1. These authors find within 
the radius range 5"^i?*^10" the J ~ H color to increase from 
1.5 to —2.3 mag, whereas the H — K color shows a continu- 
ous decrease from -0.5 mag to <-l mag. Also, the integrated 
J ~ H and H — K colors of 0.99 mag and -0.68 mag, re- 
spectively, listed in Doublier et al. differ significantly from the 
values of 0.47 mag and 0.2 mag derived in the present study. 

4.4. Henize 2-10 (ESQ 495-G21) 

This relatively metal-rich iE BCD (Zq/2A. . . '^jZq; Schaerer 
et al. '1999' Kobulnicky et al. '1999't is the first extragalactic 
system in which the broad He IIA4686 line was detected. The 
starburst nature of He 2-10 is evidenced by a chain of bright 
Super-Star Clusters (SSCs) in its brightest western SF region 
W (Conti & Vacca 1994 1, extended X-ray and Ha emission 
(Hensler et al. 119971 Papaderos & Fricke ^998 1, as well as a 
large bipolar outflow from the SF region (Papaderos & Fricke 
1998 1, with an expansion velocity between —250 and >360 
km s-^Mendez & Esteban [T9991 Johnson et al. 200^. 

Unsharp masking (inset in Fig.|6t) reveals a wealth of mor- 
phological features in the nuclear region of He 2-10, most no- 
tably an extended feature protruding northwest of region W 
(labelled NW) and an arc-like chain of compact sources con- 
necting the tip of region NW with the secondary SF knot E. 
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Fig. 4. Haro 14 (£'=12.5 Mpc). For explanations of symbols and labels, refer to Fig.|2] a): J band image and isophotes. Note 
that the center of the star-forming regions is offset by ~ 0.5 kpc from the center of the outer circular isophotes. The upper-left 
and upper-right insets show unsharp-masked versions of the central part of the BCD, marked in the main image by the white 
and black brackets, respectively. In the upper-left inset we mark regions a through d. The upper-right inset reveals an extended, 
curved chain of knots traceable up to 1 .5 kpc from the brightest region a. b),c): Surface brightness and color profiles. The thick 
grey line shows a fit to the LSB component using a modified exponential distribution Eq. (|6} with b,q ^ 3.6, 0.94 . 



Ground-based B — R maps by Papaderos & Fricke fT998^ in- 
dicate that region NW and the concatenation of sources bend- 
ing southwards of it are considerably bluer than the underlying 
LSB host galaxy. 

Our SBPs (Fig.|6]3) show an exponential intensity decrease 
in the LSB component (i?*^20"), with a J scale length iden- 
tical to that obtained previously from optical data (a = 0.67 
kpc, Papaderos & Fricke il998i) . 

4.5. Tol 1400-41 1 (NGC 5408) 

This cometary il BCD is a member of the Cen A group (van 
den Bergh 2000 1, situated at a projected distance of 7.2° (^ 
0.6 Mpc) from Cen A. Intense SF activity is taking place 
mainly at the SW part of its elongated stellar LSB host, where 
ground-based NIR and archival HSTAVFPC2 optical data (PI: 
P. Seitzer, GO-8601) reveal three bright (my ^16.2) stellar 
clusters spread over 8" (^200 pc). The appreciable ionized gas 
contribution in the SF component of Tol 1400-41 1 is reflected 
on the very blue V—I color of ^-0.4 mag, in the brightest clus- 
ter, A (cf. l.h.s. blow-up of Fig.EJi) and the large EW(H/3 )^250 
A determined for this system by Masegosa et al. J1994> . 

Interferometric 21cm line studies (Fritz l2000t show that 
the BCD is immersed within a large, rotationally supported 
H I cloud of about 8x16 kpc in size. The H i halo reveals two 



surface density maxima, one coinciding with region A and the 
other one located ^1' to the east. 

The present data, as well as Ha images by Fritz J2000> . 
reveal signatures of low-level SF activity all over the 22 J 
mag/n" isophote of the BCD, i.e 2.6 kpc across. Interestingly, 
30" NE of region A our data show a nearly circular region, 
w 12" in diameter, with relatively blue (B — R of ^0.75, 
V ~ I ^0.45 mag) colors in its unresolved interior (denoted B 
in Fig. Eh)- Such colors, being significantly bluer than those in 
the LSB host (B - i? «1, V - 7=0.6-0.7), suggest that region 
B is comparatively young. This is also suggested by the NIR 
colors (J - K=OJ. . . 0.9) and luminosities (Mjt=-8.6. . . -9.7) 
of bright sources therein which are consistent with a popula- 
tion of young (10-20 Myr) red supergiants (cf. e.g. Elias et al. 
[T9851 BerteUi et al. [T994l i. 

Another conspicuous feature seen at the northeastern part 
of the BCD is a comparatively blue (B — R ^ 0.6) curved 
strip (region C), apparently bending from the NE tip of the 
LSB component to the north; a counter-feature of this region 
is probably present at the southern part of the BCD. 

The large extent of the SF component, as well as several 
bright foreground stars in the periphery of Tol 1400-41 1, render 
the determination of the structural properties and color of its 
LSB component difficult. The SBPs (Fig.Eb) are approximated 
best with Eq. Q with (6,g)=(3.0,0.8). The mean LSB colors, 
J - H=0A2, H - Ks=0.2S and V - 1=0.6-0.7 can be brought 
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Fig.5. UM 461 (I?=14.3 Mpc). For explanations of symbols and labels, refer to Fig.|2] a): J band image and isophotes. Bright 
stellar assemblies in the central portion of the BCD are indicated. b),c): Surface brightness and color profiles. The thick grey line 
shows a fit to the LSB host galaxy using Eq. (jSJl with parameters b,q ^ 2.3, 0.85. 




Fig. 6. Henize 2-10 (D=8.7 Mpc). For explanations of symbols and labels, refer to Fig.|2l a): J band image and isophotes. The 
inset shows an unsharp-masked version of the nuclear region of the BCD (region marked by the brackets in the main image). 
b),c): Surface brightness and color profiles. 
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Fig. 7. Tol 1400-41 1 (13=4.8 Mpc). For explanations of symbols and labels, refer to Fig.|2] a): J band image and isophotes. The 
insets show magnifications of the bright SW star-forming knot A as well as the region B extending northwestwards from the 
major axis. The bluer NE region C is marked. b),c): Surface brightness and color profiles. The thick grey line shows a fit to the 
host galaxy using Eq. (|6} with b,q = 3.0, 0.82. 



into rough agreement, given the uncertainties discussed in Sect. 
13.51 The GALEV model (see Sect. 15 .4> yields for a metal-poor 
stellar population forming in a single burst or continuously an 
age between several 10^ to a few 10^ yr (cf. Fig.ll7>. However, 
these colors might be influenced by ionized gas emission and 
the younger stellar populations in the NE boundary. 

4.6. Pox 4 

Despite its moderate metal deficiency (Zq/9.3. . . Zq/7.6; 
Kunth & Sargent 1983 Vacca & Conti [T9^ . Pox 4 is a rel- 
atively luminous (Ms=-18.81) il BCG (Mendez & Esteban 
[1999) . It is accompanied by a faint SF dwarf galaxy. Pox 4B, 
at a projected distance of 5 kpc and a velocity difference of 130 
km s"'^. On our NIR images, we do not detect any emission in 
between Pox 4 and Pox 4 B down to an approximative surface 
brightness level of 24 J mag/n"(Fig.|8t). 

The optical morphology of Pox 4 is dominated by SF re- 
gions and extended ionized gas emission on a spatial scale of 
5 x2 kpc (cf. Figs. 1 and 2 in Mendez & Esteban [T999t . A sub- 
stantial color shift due to intense (EW(Ha)=1410 A) nebular 
line emission has been observed in the brightest assembly of 
SF sources M9 (inset of Fig.|8t) for which Mendez & Esteban 
J199 9> report a blue (-0.8 mag) U — B together with an ex- 
tremely red (h-0.7 mag) B — V color. The NIR colors of re- 
gion M9 (J - Hi=0.28, H - Ks^=0A7) are also suggestive of 
strong ionized gas contamination. The latter is also reflected on 
the color profiles (Fig.|8j;), approaching for R*^3" (0.7 kpc) 



colors as extreme as J — H=0.2 mag together with H — Ks=0.6 
mag. 

The SF morphology in Pox 4 is intriguing. Unsharp mask- 
ing reveals eastwards of region M9 a ring-like distribution of 
compact sources with a projected size of 3.2x1.5 kpc and a 
position angle of ^124°. A similar morphology has been re- 
ported from optical data by Mendez & Esteban (1999 1, who 
interpreted Pox 4 as a Cartwheel-like galaxy, downscaled by 
1-2 orders of magnitude. 

The available deep NIR data allow us to detect at faint 
intensities («22.5 J mag/D") a smooth underlying LSB host 
galaxy. Its Jlsb profile can be approximated with Eq. (|6j and a 
depression parameter q=0.9 (Fig.|8j3). From profile decompo- 
sition we infer the absolute J magnitude of this stellar host to - 
18.2 mag, which translates to M b >-17.2 mag for aB — J >1. 
Thus, despite its high total luminosity. Pox 4 qualifies by the 
absolute magnitude of its LSB host as a dwarf galaxy. 

Our SBPs show some similarity to the uncalibrated pro- 
files in Telles et al. (1997i out to R* = 15". However, the 
scale length derived here for 10"<i?*<16" for the LSB com- 
ponent (a=3'.'8=0.86 kpc) is twice as large as that in Telles et 
al. a997j . 

4. 7. Pox 4 B 

This intrinsically faint (Mj =-16.2 mag) galaxy has been sug- 
gested by Mendez & Esteban ( 1999 1 to have triggered the star- 
burst activity in Pox 4 through a face-on collision. As evi- 
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its companion Pox 4B 
{D=46J Mpc). For explanations of symbols 
and labels, refer to Fig. |2] a): J band image 
and isophotes. The inset shows an unsharp- 
masked close-up of the SF regions in Pox 4 
(area within brackets in the main image). M3, 
M6 and M9 mark the brightest H II regions re- 
ferred to as 3, 6 and 9 by Mendez & Esteban 
(02991. The ellipse of small boxes marks the 
ring of SF regions described by the same au- 
thors. In our data, M9 resolves into 3 knots 
arranged over a projected length of 3'.'2 (750 
pc). b)-e): Surface brightness and color pro- 
files of Pox 4 ( b),c) ) and Pox 4B ( d),e) ). 
The thick grey lines show decomposition fits 
to the stellar LSB host galaxies by means of 
modified exponential distributions (see Sect. 
K2\ . with b,q = 2.7, 0.90 for Pox 4 and 
b,q = 2.0, 0.80 for Pox 4B. 




denced by faint Ha emission (Mendez & Esteban '1999 1, Pox 
4 B still maintains a mild SF activity. The latter is probably tak- 
ing place in three compact (^Of.'S) sources, discernible in the 
central part of the galaxy. Pox 4 B shows little morphological 



distortions, except for a slight eastward extension of its LSB 
component (Fig.|S^). 

The SBPs (Fig. mi) of Pox 4B show for 2"<R*<5" an 
exponential fall-off with an a «0.29 kpc, and a flattening 
for R*<2". We have verified that the latter intensity regime 
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is not due to seeing. For this purpose, we convolved artifi- 
cial exponential 2D models with the observed a and extrapo- 
lated central surface brightness (/iE.o ~19.5 J mag/D") with 
a Gaussian kernel matching the PSF of the coadded images 
(« 0'.'7 FWHM). The resulting SBPs deviate only sUghtly from 
a pure exponential, even when the FWHM is further degraded 
by a factor of 2, and can by no means reproduce the strong 
flattening we detect in the SBPs of Pox 4B. In addition, a per- 
fect exponential distribution with the /iE,o and a quoted above 
would yield a 16% larger luminosity than the one actually mea- 
sured for Pox 4B. 

The available data do not allow us to study the NIR colors 
in spatial detail. The mean J — H and H — Kg of respectively 
RiO.4 and saO.2 are comparable to those of the LSB component 
of Pox 4. 

4.8. Tol65 

Star formation in this very metal-deficient BCD (Zq/24; Kunth 
& Sargent il983i Masegosa et al. [T994l Izotov et al. 120011 1 is 
taking place mainly in a chain of five compact (1'.'2^ ^2") 
sources located at the NE part of an irregular, blue LSB enve- 
lope (£99 1. The combination of an extraordinarily blue J — 
(~0.1 mag) with a red H — («0.8 mag) color in the sur- 
roundings of the brightest SF region e (Fig.|9t) points to a sub- 
stantial contribution of ionized gas emission. This is also sug- 
gested by the combination of a blue U — B (^-1.05 mag) with 
a moderately red B — R (^0.4-0.5 mag) color derived by P99 
all over the NE half of Tol 65, several 100 pc away from the 
opposite tips of the SF chain. Further evidence for intense ion- 
ized gas emission is provided by Keck spectroscopy by Izotov 
et al. (|2001), who derived a large EW(Ha) of > 1000 A in the 
SF component of the BCD. 

Ionized gas contamination appears to be small midway in 
the SF chain. The B-i?=0.15 mag together with J-H^ 

and H — Ks^ of 0.44 and 0.2 mag, respectively, observed in 
region C are consistent with a single burst stellar age of ^100 
Myr Note that the J — H color in region C is barely bluer than 
the average value for the LSB component (0.48 mag). This is 
also true for the H — Kg color, which could be constrained at 
R* k^Pis to <0.2 mag (cf. Fig.|9j;). 

Deep optical surface photometry (/i >28.5 B mag/D") for 
Tol 65 has first been presented in lP99l These authors found op- 
tical SBPs to show an outer exponential regime for R* >7" and 
a flattening relative to the exponential fit inwards of i?*sa3". In 
order to adequately decompose the optical SBPs, they modelled 
the LSB component with Eq. (|6} and a depression parameter 
q ^ 0.8. The type V profile of the LSB component is better 
visible in NIR wavelengths, where an even stronger flattening 
is required (q «0.9) to fit the data. The J SBP shows in the 
radius range 5'.'0 < R* < 7'.'3 an exponential intensity fall-off 
with a scale length a — 0.28 kpc, in good agreement with the 
value derived in P99 (a = 0.26 kpc). The average LSB colors 
of J - iJ=0.48±0.17, together with the a B - i? ^ 0.6 mag 
and B — J ^1.3 mag, inferable from the SBPs in lP99l are con- 
sistent with a stellar age ^ 1 Gyr, assuming an instantaneous SF 
process (cf. Sect. 13. 5> . 



4.9. Tol 1214-277 

This metal-poor (Zq/25;'F01I Izotov et al. 1200 H cometary il 
BCD undergoes strong SF activity at the NE end of an elon- 
gated stellar LSB body with an apparent size of w6.9x3.4 kpc. 
The brightest SF complex (labelled a in Fig. II Oh ) contributes 
nearly one half of the BCD's optical light within the 25 B 
mag /□" isophote jFOTl . The J* magnitude determined for re- 
gion a, 18.6 mag, shows that in the NIR the starburst still pro- 
vides ^1/3 of the galaxy's total flux. The colors of this region 
(J — H^=-Q.Q2 mag and H — Ks^=Q-5l mag) are consistent 
with a young burst age of <7 Myr for the BCD (see alsolFOU. 

Only the J band SBP of Tol 1214-277 could be observed 
with a sufficient quality for a decomposition (Fig.llOb). Similar 
to Tol 65, it shows a type V distribution, as already reported 
from deep optical VLT FORS I data by IFOTI Fitting Eq. 
to the Jlsb SBP, with the (b,q) parameters fixed to the values 
derived in FOl (3.3,0.92), we obtain a scale length a=0.53 kpc, 
slightly larger than that in the optical (a =0.49 kpc: IFOlt . 

Fricke et al. (2001) reported for Tol 1214-277 nearly con- 
stant U — B and B — R colors of respectively -0.42 mag and 
0.34 mag over an intensity span of 8 mag (colors adapted to the 
galactic absorption assumed here). Such colors, together with 
the J — H «0.4 mag and B — J «1.1 mag derived here, are 
slightly bluer that those of Tol 65, supporting the hypothesis 
that Tol 1214-277 is a relatively unevolved dwarf galaxy. 

4. 1 0. Mkn 1 78 (UGC 654 1 ) 

The spectroscopic properties of this intrinsically faint (Mb=- 
13.9 mag) iE BCD, a member of the CVn cloud I of galax- 
ies (Makarova et al. I1998> . have been investigated in e.g. 
Gonzalez-Riestra et al. (1988l and Guseva et al. (2000 1. Optical 
images reveal a complex morphology in the SF component, 
notably a pronounced separation of HSB regions into a com- 
pact assembly of bright knots to the SE (regions a and b) and 
an extended arc-like segregation of fainter sources to the NE. 
NIR images, on the contrary, show a more regular morphology, 
with only a few prominent sources in the central part of Mkn 
178 (see Figll IhV Region a (toj-'-=17.5 mag; cf. inset in Fig. 
Illh ) coincides with the optically brightest region in the BCD 
and is the main locus of active star formation. The brightest 
NIR source, b (m,7-'-=16.5 mag), is optically faint and shows 
nearly the same B — R color as the surrounding diffuse emis- 
sion within the plateau component (cf. Papaderos et al. l2002t . 
Sources C and d are immersed within the northern featureless 
region. The non-detection of an optical counterpart for region 
(j, and the overall optical/NIR morphology of the SF compo- 
nent are suggestive of inhomogeneous, large-scale dust absorp- 
tion on a spatial scale of '-^l kpc. This may cause the apparent 
separation of the SF component in two large detached com- 
plexes and hide knot d in optical wavelengths. 

Gonzalez-Riestra et al. (I1988> and Guseva et al. (I2000> 
have shown ionized gas emission to be negligible in region b 
(EW(H/3)=24. . . 34 A). By the col or exc ess of E(B - y)=0.25 
given in Gonzalez-Riestra et al. ( 11988! ). its observed B — R 
color (0.67 mag, Papaderos et al.'2002'l transforms to 0.28 mag, 
suggesting a single-burst age of ^100 Myr. This is consistent 
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Fig. 9. Tol 65 (£'=34.2 Mpc). For explanations of symbols and labels, refer to Fig.|2] a): J band image and isophotes. The sources 
b, C and e along the NE-chain of SF regions are marked, following the notation of P99 The feature marked by an asterisk is 
a red (J — Ks ~ 0.8 mag) compact source probably not belonging to the galaxy. The faint (mj « 20.8 mag) galaxy "G1 " (cf. 
IP99 ') is visible to the SW. b),c): Surface brightness and color profiles. The thick grey line shows a decomposition fit to the host 
galaxy by means of a modified exponential distribution (see Sect. l3.2l i with b,q = 2.7, 0.90 . 
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Fig. 10. Tol 1214-277 (L»=102.6 Mpc). For explanations of symbols 
and labels, refer to Fig.|2l a): J band image and isophotes. The bright 
star-forming complex at the NE edge of the galaxy (a) is marked, along 
with the faint galaxies G1-G3 (see FOlj . b),c): Surface brightness and 
color profiles. The thick grey line shows a decomposition fit to the host 
galaxy by means of Eq. ^ with 5, g = 3.3, 0.92 . 




Noeske et al.; NIR observations of Blue Compact Dwarf Galaxies 



19 




4U U -iiO -4G -UO 10 10 2o 3o 4o so 

arCO-eC: R [arcsec] 



Fig. 11. Mkn 178 (D=4.2 Mpc). For explanations of symbols and labels, refer to Fig.|2] a): J band image and isophotes. The 
inset shows a close-up of the star-forming regions (delimited by brackets in the main image). The regions a-d described by 
Gonzalez-Riestra et al. ( I1988I I and Papaderos et al. J2002> are labeled. b),c): Surface brightness and color profiles. 



with the colors of J - * =0.24, H - A"* =0.06 and B - J* = 1 .6 
mag obtained in the present work, if they are de-reddened 
adopting the same amount of intrinsic extinction. The presence 
of stellar complexes with an age of the order of '--^lO^yr is in 
line with the detection of numerous luminous AGB stars in NIR 
color-magnitude diagrams (Schulte-Ladbeck et al. 120001) . from 
which these authors infer significant SF activity over the last 
few 10^ yr. 

NIR SBPs of Mkn 17 8 (Fig.[TTt )1 resemble closely the op- 
tical ones (Papaderos et al."2002^, showing an extended plateau 
component ontop a smooth, exponential LSB host galaxy. For 
the latter we derive a J band scale length of 0.28 kpc and colors 
(J - H «0.53, H - K' «0.2, B - J «2.0 wA B - R «1.1 
mag) that consistently indicate an old stellar^ LSB background. 

4. 1 1. Mkn 1329 (IC 3589/91, UGC 7790, VCC 
1699) 

This relatively metal-rich {Zq/5, Guseva et al.'2000') cometary 
il BCD (Mb ~ -16.8, Yasuda et al. .l997,) is associated with a 
group of 1 1 galaxies (LGG 296, Garcia [T993> within the Virgo 
Cluster. 

As shown in the contrast-enhanced image (inset in Fig. 
I12h). numerous irregular concentrations are present in the inner 
zone of its elongated stellar LSB component, over a projected 
length of ~3.5 kpc. An Ha exposure by Gallagher & Hunter 
J1989> shows that SF activity is almost exclusively confined to 
a bright H II region (marked with a cross in Fig.ll2h') at the SW 
end of the high-surface brightness component. The Ha flux de- 



rived for that SF region by the latter authors yields, at the dis- 
tance adopted here, an Ha luminosity of «1.5xlO^" ergs^^ 
and a star formation rate (SFR) of '--^ 0.1 Mq yr~^. The dom- 
inant SF region in Mkn 1329 is identifiable with a compact 
(© ^ 300 pc), luminous (A/j = -14) NIR source; its colors, 
(J - H)i=0.l6 and (H - Ks)^=0.5l, ai-e probably significantly 
influenced by nebular emission (EW(H/3)=276 A, Guseva et al. 
20001. 

The morphology of the SF regions and the high luminosity 
of the SW complex in Mkn 1329 are typical among cometary il 
BCDs. The similarity of these systems with respect to the mor- 
phology of their LSB and SF component suggests a comparable 
evolutionary state/history or common mode of star formation. 
Propagation of SF activities along those objects' major axes 
has been proposed and observationally supported by a number 
of recent studies (cf. Noeske et al. '2000 and references therein). 

The SBPs of Mkn 1329 (Fig. 112b') indicate a moderate cen- 
tral flattening of the LSB component, which can be approxi- 
mated both by a Sersic law with 77lsb ~ 0.7, or a med with 6, q 
= (1.6,0.7) (see Tableland Sect. 13. 2> . The detection of a mod- 
erate central flattening in the LSB component owes to the com- 
paratively small contribution of the superposed younger stel- 
lar population in Mkn 1329 (wl7% of the total J band light). 
The NIR colors in the outskirts of Mkn 1329, (J - i7)*=0.64, 
(H — A's)*=0.06, suggest an evolved stellar LSB component. 
Adopting a B — J color characteristic of a few Gyr old stel- 
lar population, we can estimate the B band structural param- 
eters from the /iE,o, « and Mlsb determined from J SBPs. 
We find that the host galaxy of Mkn 1329 shows, similar to 
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Fig. 12. Mkn 1329 {D=16 Mpc). For explanations of symbols and labels, refer to Fig. |2] a): J band image and isophotes. 
inset: /ife-transformed image of the central area (delimited by the brackets in the main image), showing the distribution of 
the bright regions along the major axis. The white cross marks the bright SW star-forming region for orientation. b),c): Surface 
brightness and color profiles. The thick grey line shows a decomposition fit to the host galaxy by means of a modified exponential 
distribution, Eq. (|6}, with b,q — 1.6, 0.70 . 



other cometary BCDs (Noeske et al. I2000> . structural proper- 
ties intermediate between extended dIs/dEs and compact iE/nE 
BCDs. 



4.12. IC4662 (ESO 102-G014, He 2-269) 

This nearby {D=2 Mpc; Heydari-Malayeri et al. '1990) dwarf 
galaxy has been studied spectroscopically by e.g. Pastoriza & 
Dottori J1981t . Stasinska et al. J1986> . Heydari-Malayeri et 
al. (199^) and Hidalgo-Gamez et al. ('200ri. The latter au- 
thors derived the oxygen abundance in the brighter SF re- 
gions A&B to Zq/6.5. . . Zq/7.6, respectively. A continuum- 
subtracted Ha map (Fig. 113b) from Papaderos et al. (2002) re- 
veals an extended and complex morphology of the ionized gas 
emission in the upper half of the galaxy, where the Ha emis- 
sion peaks, and shows a number of shells extending up to '^0.5 
kpc NE. Another interesting feature is the extranuclear region 
D, seen in the outer regions (/i ^ 22 J mag/D") of IC 4662. 
Hidalgo-Gamez et al. ("2001 ) find this Ha emitting region to be 
less metal-rich than the central SF regions (A), and to show 
a significant recession velocity difference (250 ±150 km/s, 
Hidalgo-Gamez et al. 2002 ) to IC 4662. Our deep images show 
that region D is not well detached from the main body of dif- 
fuse Ha emission in IC4662 (contrary to Heydari-Malayeri et 
al.|l,990), but apparently connecting with regions A&B through 
a chain of Ha sources. 



Hidalgo-Gamez et al. J2002> have proposed that region D 
may be either a chemically and kinematically distinct complex 
within IC 4662, or a close companion object. The latter possi- 
bility is particularly interesting in view of the hypothesis that 
very close, gas-rich dwarf companions might be conceivable 
triggering agents of starburst activity in BCDs (Taylor et al. 
rT995l Pustilnik et al. l200Tl Noeske et al.|200D- 

From the available data, we resolve a wealth of morpholog- 
ical information in the central part (Fig.ll3h) of IC 4662, no- 
tably the head-tail morphology of the starburst. Luminosities 
and colors of the brightest point sources in region A are in 
agreement with those of red supergiants, supporting the results 
by Heydari-Malayeri et al. J1990t . as well as of blue super- 
giants and AGB stars (cf. the upper left inset in Fig.ll3h: see 
Schulte-Ladbeck et al. . 2001 for color limits separating the lat- 
ter classes of giant stars). 

As IC 4662 is located at low galactic latitude (-17.8°), pho- 
tometric studies of its LSB component are complicated by the 
dense foreground Galactic stellar field. Also residuals in the 
removal of the NW bright star may affect the photometry for 
faint isophotal levels. A tentative exponential fit to the J SBP 
yields for i?*>65" , i.e. outside significant nebular emission, 
a scale length of ~0.21 kpc. However, inspection of the Jlsb 
profile shows that the BCD follows a type V distribution, fitted 
best with Eq. (|6} with a central surface brightness of 16.8 J 
mag/n", a scale length of ~150 pc and a depression parameter 
as large as q fa 0.97. If so, the starburst contributes ~50% of 



Noeske et al.; NIR observations of Blue Compact Dwarf Galaxies 



21 



KJO 



")0 - 



'.0 



at, ^ 



-100 



IC«52 
J band 



(a) 




..'.1- 





• 










DXTipancrrt D 



■ I I 



■l;o 







100 





Fig. 13. IC 4662 {D=2 Mpc). For explanations of 
symbols and labels, refer to Fig. |2l (a): J band 
image and isophotes. The insets show magnified 
views of the star-forming region A in the cen- 
tral part of the BCD and of the D component, 
about 1 'A from A. Some bright sources in region 
A are classified according to Schulte-Ladbeck et 
al. 



the total J light of the BCD. IC 4662 was hitherto classified as 
a dwarf irregular galaxy (cf. e.g. de Vaucouleurs et al. .1991 ) . 
However, the structural properties of its LSB host, as well as 
its intense and spatially extended SF activity, place IC 4662 in 
the range of BCDs, making it probably one of the closest BCDs 
known. 



4.13. UM488 (Mkn 1304, SBS 1139+006, UGC 
6665) 

This is a distant (D « 76 Mpc, Mirabel & Sanders [T988 1 
blue compact galaxy, known to be relatively metal-deficient 
(Zq/S.S; Masegosa et al. [T994l Izotov & Thuan .l999j . Despite 



Boon into red supergiants ( J - H ^0.15, 
crosses), blue supergiants (J — H <0.5, squares), 
and AGB stars (J - H >0.9, lozenges), (b): 
Continuum-subtracted Ha map (from Papaderos 
et al. I2003t . showing that all regions (A-D from 
Heydari-Malayeri 1990 1 are Ha emitting. (c),(d): 
Surface brightness and color profiles. The thick 
grey line shows a fit to the host galaxy using a 
modified exponential distribution, Eq. (|6j, with 
6,(7 = 4.8,0.97. 
its resemblance (Fig.ll4h) to cometary il BCDs, UM 448 does 
not qualify as a dwarf due to its intrinsic luminosity {Mb ~- 
19.7 mag) and large linear extent (major axis length ^25 kpc). 
Active star formation is taking place within an extended (diam- 
eter ~7 kpc) high surface brightness region at the NE part of 
the galaxy, showing some fainter SW extension. The starburst 
component contributes ?a70% of the J band light. Guseva et 



al. J2000I I estimated the starburst luminosity to be powered by 
1.6 X 10^ O stars (value referring to the distance assumed here), 
0.6% of which are undergoing their Wolf-Rayet phase. Studies 
by Sage et al. (I1992> indicate an ongoing SFR of 21 il/© yr~^, 
and a large reservoir of molecular gas amounting to 2.8 x 10^ 
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Fig. 14. UM 448 {D=l 6.1 Mpc). For explanations of symbols and labels, refer to Fig.|2] a): J band image and isophotes. The 
inset shows a contrast-enhanced close-up of the star-forming regions (delimited by brackets in the main image); individual knots 
(a— C) are labeled. b),c): Surface brightness and color profiles. 



Mq or ~l/2 of the H I mass (5.8 x 10^ A/©; Mirabel & Sanders 
11988 1 of UM448. 

On unsharp-masked images, the SF region splits into 3 
high-surface brightness entities (see the inset in Fig.ll4h). the 
brightest one being the reddest (J — =0.7 mag) whereas the 
fainter regions b and C are slightly bluer (J — H^k,Q.5 mag). 
The SBPs (Fig. 114b) show an exponential decay for i?*>15" 
with a J band scale length of 3.8 kpc, somewhat smaller than 
the value ~5.2 kpc (14") infeiTed in Telles et al. ( fT997l . 

5. Discussion 

5. 1. The shape of NIR surface brightness profiles 

The SBPs derived for the sample BCDs in Sect.|4]bear close re- 
semblance to those typically inferred from optical broad-band 
data (cf., e.g.. |P96al Telles et al. [T997l Marlowe et al. 1997 
Doublier et al l 19971 19991 Cairos et al. 12001 at . In most cases, 
SBPs show in their outer part an exponential intensity decrease 
and red, nearly constant colors. This outermost SBP part is at- 
tributable to the underlying LSB host galaxy which, except for 
a few rare cases, shows clear evidence for a several Gyr old 
stellar population. At intermediate and small radii the emission 
of a BCD is dominated by the younger stellar component. In 
the main class of iE BCDs (LT86') its luminosity output is re- 
flected on two conspicuous SBP components. 

1) At small radii (i?*^100 pc), a feature commonly ob- 
served in optical and NIR SBPs is a central intensity excess. 
This is seen in e.g., the J SBP of Tol 3, Mkn 178, Mkn 1329, 



Haro 14 (our sample) or Mkn 36 and UM 462 (03^. This com- 
ponent is due to the brightest and typically youngest stellar as- 
sembly. In most cases, this narrow innermost excess is due to 
one and the same region in both optical and NIR wavelengths. 
An exception is presented by Mkn 178, where the brightest re- 
gion in the NIR is offset from that in the optical by ^^150 pc, 
or by the BCD Mkn 35 in the sample of lC03l 

2) The second SBP feature has been referred to as plateau 
in IP96al It shows a nearly constant or slowly decaying inten- 
sity, and extends typically out to /i ^^24 B mag/D". SBPs with 
a prominent plateau on top a more extended exponential LSB 
envelope can in general not be adequately fit by a simple func- 
tion (e.g. a Sersic profile) over their whole intensity span. A 
meaningful decomposition of such SBPs is only possible when 
an extra component (e.g., a Sersic distribution with rj < 1) is 
introduced in order to fit the plateau ( P96a' Cairos et al.'2001a1. 
There is, so far, no observational support for the plateau being 
a dynamically distinct or even interaction-induced stellar entity 
with nearly constant M/C. Quite contrary, appreciable color 
gradients (<,1.5 B-Rmag kpc"^) in the HSB part of BCDs in- 
dicate that the plateau light is mainly due to a young and mod- 
erately evolved stellar population. 'P96b have shown that a con- 
spicuous plateau in the SBPs of many iE BCDs can naturally 
result from the superposition of diffuse and compact SF sources 
with varying luminosity and galactocentric distance on a more 
extended, exponential LSB component. Conversely, the plateau 
is nearly absent in the intrinsically brighter nE BCDs, where 
most of the starburst light originates from the nuclear region of 
the BCD. The overall SBPs of these systems show frequently 
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a concave shape, fitted satisfactorily by Eq. (|8) with a Sersic 
exponent rj >1, and sometimes resemble a de Vaucouleurs pro- 
file. 

Discerning the formation history of the stellar populations 
memorized in the plateau light poses a challenge for surface 
photometry and CMD studies, as nebular line emission (e.g. Tol 
65; Papaderos et al. P99 Izotov et al. 2001 Tol 1214-277;|F0D 
and patchy dust absorption (Tol 3, Mkn 178; this paper, Mkn 
33 and Mkn 35;|C03 I Zw 18; Cannon et al. 120021 may both 
hamper standard age-dating techniques. In addition, an impor- 
tant and mostly overlooked source of systematic uncertainties 
in the determination of colors or EWs within the plateau stems 
from the unknown line-of-sight contribution of the underlying 
old LSB component. As pointed out in Sect. 13. 31 the latter may 
redden colors of compact stellar clusters in the SF region by 
up to a few tenths of magnitude. Corrections of this order have 
also been reported for the iE BCDs Mkn 370 and Mkn 178 by 
Cairos et al. (20()2a) and Papaderos et al. (2002i, respectively. 
Evidently, for such a correction one has to assume a model for 
the intensity distribution of the LSB component just beneath 
the SF regions (i.e., for i?*<f ^). It is a common practice to ex- 
trapolate the exponential slope of the LSB periphery of BCDs 
all the way to i?*=0". However, the universal validity of this 
procedure for BCDs has been questioned in, e.g., P96a. These 
authors discussed observational evidence for type V profiles in 
BCDs and proposed the alternative fitting formula Eq. ^ to 



approximate such convex distributions. Alternativelv. lC03l dis- 
cuss cases where the LSB host galaxy shows a concave profile 
which is better fit by a Sersic model with an exponent 77 > 1. 
The photometric structure of the LSB host of BCDs will be dis- 
cussed in the light of the present NIR data in the next section. 

5.2. Hints to a centrally flattening exponential 
distribution in the underlying stellar LSB 
component? 

Nine of the BCDs included in our sample (Table|2j show signa- 
tures of a type V profile in their LSB hosts. This is because an 
exponential law, or any Sersic distribution with ?/ > 1, fitted to 
the outer part of their J SBPs, predicts at intermediate to small 
radii a higher intensity than the observed value (cf. Sect. 13. 2> . 
Of course, pure exponential fits to Jlsb profiles cannot always 
be definitely ruled out within the la uncertainties. However, 
such fits would either overestimate the central intensity of the 
stellar LSB host, hence underestimate the luminosity fraction 
of the SF component (cf. Sect, \5.2i or systematically overes- 
timate the J band exponential scale length, thus imply an im- 
plausibly large B — J color gradient for the old underlying LSB 
population. 

The limited size of our present sample does hardly allow 
to estimate the frequency of type V profiles in BCDs. This is- 
sue will be addressed in a forthcoming paper of this series, fo- 
cussing on the complete NIR sample. However, the evidence 
gathered so far (see also Sect. 13. 2t strongly supports the idea 
that a substantial fraction of BCDs shows type V profiles in 
its LSB host galaxy. The rare detection of type V profiles in 
BCDs in previous optical studies can be explained by the fact 
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Fig. 15. Decomposition of the B band SBP of Haro 14 (filled 
circles) assuming for its LSB host a pure exponential model 
(dotted line) or a modified exponential distribution (Eq.|6j thick 
gray curve) with the flattening parameters (fe,q)=(3.6,0.9) de- 
rived in the J band. The surface brightness distribution of the 
star-forming component, as obtained by subtraction of the re- 
spective LSB model from the B SBP is illustrated by the filled 
and open interconnected symbols (see discussion in the text). 



that in those wavelengths extended starburst emission over- 
shines the LSB host within typically its inner ~2 exponential 
scale lengths (P96F Noeske 1999", Cairos "2000' Cairos et al. 
|2p01ai. As a result, a type V profile in iE/nE BCDs could be 
detected only when a significant depression of the exponential 
LSB slope occurs for R*> 2a, i.e. if 6 > 2 in Eq. (jSJl. NIR ob- 
servations do not overcome, but alleviate the problem of severe 
light pollution by the burst, proving in many cases crucial for 
disentangling a pure exponential from a type V distribution in 
the underlying LSB galaxy. 

This is illustrated in Fig.^jon the example of the iE BCD 
Haro 14. Its B SBP is roughly exponential out to i?*=40" 
(/I ^27 B mag/n"), suggesting that the exponential slope seen 
in the LSB periphery continues all the way to i?*=0". In that 
case, one is tempted to conclude that the ongoing burst gives 
rise to a modest luminosity increase, only. Indeed, subtraction 
of the exponential LSB model from the SBP suggests that SF 
sources (small interconnected squares in Fig. I15> are all con- 
fined within i?*=8" and that they account for no more than 7% 
(>16.5 mag) of the B light of the BCD. Evidently, this con- 
clusion is hardly compatible with the extended morphology of 
the SF component and copious Ha emission of Haro 14 (cf. 
Marlowe et al. [T997l Doublier et al. [T999l . 

One arrives at a quite opposite conclusion on inspection of 
the NIR profiles. Figure^jshows that a pure exponential law 
falls short of fitting the J SBP inwards of i?*<20"; a plausible 
fit to the NIR data must invoke a significant flattening of the ex- 
ponential slope inside ^3 scale lengths of the LSB component. 
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An adequate fit to the Jlsb is best achieved with a med (Eq.|6j 
with (b,q)=3.6,0.94. Such a strongly flattened exponential LSB 
model yields a stellar mass by a factor ^^3 smaller than that 
of a purely exponential profile and implies that the starburst 
(small interconnected circles in Fig. I15> contributes ^60% of 
the B light within the Holmberg radius, i.e. more than a factor 
of eight larger than the value inferred before. The isophotal B 
magnitude within 25 B mag/D" is then increased by ^1.1 mag 
due to the starburst. This value is still compatible to the mean 
value of 0.75 mag deduced for BCDs in P96b. and Salzer & 
Norton |l999l. 

The considerations above show that the choice of the model 
for the LSB emission may have far-reaching implications to our 
view about the SF amplitude and photometric fading of BCDs. 
The usual assumption of a purely exponential LSB model may, 
in some cases, strongly underestimate the luminosity fraction 
of the young stellar population, and lead to the conclusion that 
a "burst" is merely a minor luminosity enhancement in the life- 
time of a BCD. If, on the contrary, the LSB emission is assumed 
to follow a type V profile, then the estimated luminosity of the 
superimposed young stellar population may increase by more 
than a magnitude. 

Consequently, different LSB models imply different 
amounts of photometric fading of a BCD once SF activities 
have terminated. This has to be borne in mind when discussing 
evolutionary links between BCDs and other dwarf galaxies, as 
well as a possible fading of luminous blue compact galaxies 
at medium redshift to local spheroidals (cf. e.g. Guzman et al. 
Il998 l. Note that the assumption of a type V profile has no ef- 
fect on the extrapolated central surface brightness /iE,o and the 
exponential scale length a of the LSB host galaxy of a BCD. 
As a result it is not expected to significantly change the system- 
atic difference between BCDs and other types of dwarf galaxies 
on the /iE.o--^^LSB and log(Q;)-AfLSB parameter space I P96bl 
Marlowe et al. 1997 Salzer & Norton 1 999. Papaderos et al. 
I2OO2I . 

A dedicated investigation of the frequency and origin of 
systematic deviations from the exponential law (e.g. a type V 
profile, or Sersic profile with 77 > 1) in the stellar LSB host of 
BCDs is apparently of great interest. A first step towards such 
studies is to test the suitability of different empirical functions 
in approximating the observed LSB intensity distribution. This 
issue will be briefly discussed in the following section. 

5.3. The Sersic law vs the med for studies of type V 
SBPs 

As noted in Sect. 13.21 also the Sersic law, commonly applied 
to structurally analyze various extragalactic systems, can ap- 
proximate type V SBPs. In this section, we further discuss why 
in this work preference was given to the med model to fit and 
quantitatively study such SBPs. 

Both the Sersic law and the med are fitting functions of 
empirical origin (Sect. 13 .2> . The choice of either model does 
therefore not imply any assumptions on the physical back- 
ground of type V SBPs, such as the dynamics or stellar mass 
distribution of a galaxy. However, although either function was 
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Fig. 16. upper panel: Fit to the LSB component of Haro 
14 (thick gray curve, labelled 4) assuming an inwards flat- 
tened exponential distribution (med, Eq. |6j with the extrap- 
olated central surface brightness /iE,o and exponential scale 
length a indicated at the upper-right, and flattening parame- 
ters (&,g)=(3.6,0.94). Distributions 2 and 5 correspond to a med 
with equal /iE,o and a but different (&,g), and 1 describes a pure 
exponential distribution (q = 0). The dashed curve shows a fit 
to distribution 4 using a Sersic fitting law (Eq.|8j. lower panel: 
Residuals between the assumed med distributions 2 through 4 
(upper panel) and Sersic fits. 

found to satisfactorily approximate the type V SBPs found in 
this work, they principally differ with respect to their shape. 
This difference is illustrated in the comparison shown in Fig. 
[T6l where we fit for illustrative purposes a Sersic profile (thin 
dashed line) to the med model (heavy Une, labelled 4) obtained 
for the LSB component of Haro 14 in Sect. |3 Distributions 
2 and 3 illustrate med profiles with equal extrapolated cen- 
tral surface brightness /iE,o (17.46 J mag/D") and exponential 
scale length a (0.37 kpc) as 4, but differing degrees of flatten- 
ing (q) with respect to the pure exponential (q — 0) model 1. 

The med 4 compares best to a Sersic profile with a cen- 
tral surface brightness /is,o = 20.45 J mag/D", a scale length 
(3 «1.14 kpc and an exponent rj — 0.57. Note that the deduced 
fis.a compares well to the actual central surface brightness of 
the LSB host galaxy 

/iE,o + 2.5 log ~ "^^-^ mag/D", (9) 

implied by the med model. 

It is evident that the Sersic law and the med distribution in 
Fig-Elclosely follow each other over a radius range of nearly 
10 scale lengths a. Nevertheless, the lower panel of Fig. ^| 
shows that the residuals between either model are small (^0.2 
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mag), but systematic. This difference reflects the fact that a 
Sersic law with 77 < 1 differs from an exponential at all radii, 
while the med approaches an exponential slope for R > ba. 

We recall that both a Sersic and a med model were found 
to approximate well the inner portion of a med profile (e.g. 
within ^1.5 cutoff -radii ba). The question which distribution 
gives a more adequate description of the observed type V SBPs 
can hence only be assessed through studies of the outer parts 
of such profiles. If the type V SBPs show a perfectly exponen- 
tial fall-off in their outer parts, continuing several scale lengths 
beyond a cutoff radius ba (cf. Eq. |6}, a Sersic law will fail 
to correctly reproduce these light distributions, and no Sersic 
model fit to the LSB emission can correctly recover the expo- 
nential scale length a and will not be free of ambiguity (Sect. 
15.3. 1> . Alternatively, if subsequent studies show that the LSB 
profiles deviate systematically from the exponential law at all 
radii, then a Sersic profile might be the preferable fitting law, 
albeit the drawbacks discussed in the following Sect. 15.3.11 

Such considerations call for deeper photometric studies of 
the LSB component of BCDs, and in dwarf galaxies in general. 
An ultimate check requires deep surface photometry, allowing 
to derive SBPs out to ^ 10a with an accuracy better than 0.2 
mag (see the residuals in Fig. I16> . not reached by currently 
available data. 

As long as tests of this kind await to be done, we argue, in 
agreement with XIQ3 that a pure (Eq.|5} or modified exponen- 
tial (Eq. |6j fitting formula should be preferably used to fit the 
underlying LSB component of BCDs. Previous optical studies 
(e.g.jJ M. P96a Telles et al. [T997l Marlowe et al..l997. Salzer 
& Norton [T999l Vennik et al. l200()l Cairos et al. l20Qlall2()Qlbl 
Makarova et al. "2002'l of BCDs suggest that those systems 
mostly show no strong deviations from an exponential law over 
their outer LSB SBPs. Also if small deviations from an expo- 
nential law which are hardly detectable in current data should 
exist, a med gives a robust approximation by fixing the LSB 
profile shape. While a Sersic law at first view offers a greater 
flexibility in this profile region, its application in the extremely 
low S/N regime to quantify marginal deviations from an ex- 
ponential is problematic, as discussed in the following section. 

5.3.1 . A cautionary note on structural analyses of 
BCDs using the Sersic law 

A drawback of Sersic models when fitted to type V profiles is 
that the solution can vary considerably, depending on the fitted 
radius range. For instance, a Sersic fit to the med 4 (Fig.ll6> 
for <23.5 J mag/D", or within 5a (--30") yields an 77 = 0.45, 
significantly lower than the value above. Note a Sersic expo- 
nent rj <0.5 reflects a central minimum in the radial luminosity 
density distribution of the LSB component (P96a). The stabil- 
ity of the shape parameter 77 has been explored in detail in C03l 
these authors remarked that a Sersic fit can be very uncertain 
for concave SBPs, where the exponent 77 can vary by up to an 
order of magnitude, depending on the fitted range and the radial 
sampling of points in the profile. 

Another objection to the view that Sersic models offer a 
robust tool to systematize the structural properties of dwarfs 



comes from the degeneracy of the shape parameter rj vs. the 
pseudo-scale length f3. A fit of Eq. (|8j to the med distribution 
2 yields an 77=0. 94 and a (3 «0.44 kpc (> a). As for profile 3, 
it is best approximated with an 77 = 0.75 and (3 « 0.72 kpc 
i2a). We see that while the Sersic exponent 77 decreases with 
increasing degree of flattening, i.e. as we go from distribution 2 
towards 4, the corresponding scale length [3 increases from 0.44 
kpc to 1.15 kpc. Actually, the exponential scale length a of the 
med profiles considered here is nowhere recovered by fitting a 
Sersic model. Instead, one is left with the pseudo-scale length 
(3, which by itself alone carries no quantitative information on 
the structural properties of the LSB component. 

Obviously, the scale length (3 is only meaningful in con- 
nection with rj. However, these two quantities are rendered 
impractical for a systematic study of the LSB component by 
their strong non-linear coupling (see discussion in e.g. Young 
& Currie 1994 Cellone & Buzzoni 2001 1, and dependence on 
the fitting procedure (see above). Furthermore, observational 
uncertainties connected with, e.g., the sky subtraction (see 
e.g. Ceflone & Buzzoni 200l1 IC03t . filtering of images prior 
to surface photometry, and imperfect removal of background 
sources may also skew the f3 vs. rj parameters in a hardly pre- 
dictable manner, making Sersic fits to the outer part of SBPs a 
hazardous procedure JC03t . These problems are probably not 
worrisome in studies of early type galaxies, i.e. systems with 
little morphological distortions and a nearly constant A4/C, 
where also the central, high S/N regions are accessible to con- 
strain the global photometric structure. In irregular SF galax- 
ies, however, the derived Sersic parameters depend on sub- 
jective choices in the data processing (see above), profile fit- 
ting and, equally important, on the profile extraction methods 
themselves. This is particularly important in the case of BCDs, 
where mostly only the outer, low S/N part of the LSB host is 
accessible to structural studies. 

It is conceivable that the latter circumstances imprint the 
Sersic parameters ifJ.s,o,P,v) inferred recently for a sample 
of luminous BCGs by Bergvall & Ostlin (l2002l . We be- 
lieve that the extremely large 77LSB exponents they derive 
(up to 20), in conjunction with central surface brightnesses 
<C Omag/arcsec^, should be regarded formal only, rather than 
a manifestation of an extraordinarily dense Dark Matter halo 
which, as advocated by these authors, dominates the stellar dy- 
namics of a BCD. 

5.4. The colors of the stellar LSB host galaxy 

The available NIR data allowed not only for the detection and 
structural analysis of the LSB component for all sample BCDs, 
but also for the determination of at least one NIR color in the 
LSB host (Sect. |33}. 

The B - J vs. J - H and H - Ks vs. J - H LSB col- 
ors (Table|3}, whenever available, are displayed in Fig.[n] For 
comparison, color predictions from the GALEV evolutionary 
synthesis model (Schulz et al. 2002, Anders, Bicker & Fritze - 
V. Alvensleben 2002, priv. comm.) are overlaid in both panels, 
calculated for stellar populations with metallicities of Zq/50 
and Zq/2.5, formed either in an instantaneous burst or continu- 
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Fig. 17. Colors of the stellar LSB host galaxies (filled circles) as listed in Table |3l For comparison, we show the color evolution 
of synthetic stellar populations, calculated with the GALEV evolutionary synthesis model (Schulz et al. 2002 Anders, Bicker & 
Fritze-v.Alvensleben 2002, priv. comm.). An initial mass function with a Salpeter slope and respective lower and upper stellar 
mass cutoffs of 0.08 and 100 M0 was assumed. We show single-age stellar populations, formed in an instantaneous burst, as well 
as in a continuous process with a constant SFR, both calculated for metallicities of ZqI5Q and Zq/2.5. For clarity, only population 
ages of 0.1, 0.5, 1, 3, 7, 10 and 14 Gyr are shown. The youngest (0.1 Gyr) point of each evolutionary track coiTesponds to the 
bluest H — Ks and B — J color. All NIR colors shown are transformed to the 2MASS photometric system (see Sect. 12.4.11 . while 
the B band refers to the Johnson system. Note that the eiTor bars indicate cumulative uncertainties due to all major eiTor sources 
(see Sect. 13. 5> . The colors of ionized gas emission from a metal-deficient (Zq/20) H II region (Te = Iff* K, He =100 cm , 
Case B recombination; cf. Kruger[l992j are indicated by the asterisk. Arrows show reddening vectors. 



ously with a constant SFR. As evident from Fig.^] NIR colors 
(J — H ws. H — Ks) show little evolution for old stellar popula- 
tions (^0.1 mag for an age >0.5 Gyr). However, optical-NIR 
colors, such as, e.g., B — J (right panel of Fig. ll7> . change by 
up to ~1.5 mag for ages >0.5 Gyr. Within the above discussed 
errors, such colors will allow to constrain the age of the stel- 
lar LSB host of old lE/nE BCDs with a precision of 1-2 Gyr, 
once precise spectroscopic measurements of its metallicity be- 
come available (see the discussion of Gil de Paz et al. .2000b) . 
Fig- El shows that the observed colors are generally compati- 
ble with the GALEV predictions within their Icr uncertainties. 
This is also the case for model predictions from PEGASE (Fioc 
& Rocca-Volmerange[l997 1. 

Five BCDs in our sample, Pox4, Pox4B, Tol 1400-411, 
Tol 65 and Tol 1214-277 (TableQ, show blue ( J ~ H <0.5, 
B — J <1.5) colors in their LSB component. This may be 
partly due to an appreciable contribution of extended nebular 
emission, as suggested by previous optical spectrophotomet- 
ric work. The blue optical-NIR LSB colors of Tol 65 and Tol 
1214-277 maybe attributed to a comparatively young (^1 Gyr) 
photometrically dominant stellar population (cf. Papaderos et 
al P99 , F01 1. The optical-NIR colors of the LSB host of the 
remaining sample BCDs can be reconciled with an old stel- 
lar population of subsolar metallicity. This result is in line 
with most published optical-NIR photometry on stellar hosts 
of BCDs (Gil de Paz et al. l2000bl Vanzi et al. l2002l Papaderos 
et al. '2002), which indicates a comparable range of colors for 
the LSB population of BCDs. 



We can compare the NIR colors of the LSB component 
with data from the literature for three sample BCDs, only (Tol 
3, Haro 14 and UM 461). The colors of Tol 3 compare well 
with those reported by Vanzi et al. (I2002> over the whole ra- 
dius range of (R*=0". . .50 ")• Our photometry for Tol 3 and 
Haro 14 is also in agreement with Doublier et al. (2 001i) for 
small radii (for i?*^15" and i?*<10", respectively). For larger 
radii, however, the NIR colors by Doublier et al. are not com- 
patible with our results within the la uncertainties. Outside the 
star-forming component the color profiles by Doublier et al. ap- 
proach values of J — H > 2 for Tol 3, and J — i7 w 1.0 . . . > 
2.5, - « ... - 0.6 for Haro 14. For UM 461, their color 
profiles vary between « 1.1 ... 2.3 (J — H) and « —0.6 ... 0.5 
(H - K). 

6. Summary and Conclusions 

We have analyzed deep Near-Infrared (NIR) J, H, K images 
of 12 Blue Compact Dwarf (BCD) Galaxies and one lumi- 
nous Blue Compact Galaxy. These objects, together with those 
studied in an accompanying paper (Cairos et al. 2003), con- 
stitute the first part of a sample of 40 BCDs for which deep 
NIR images were obtained in the framework of a large-scale 
multi-wavelength study. The limiting surface brightnesses of 
our data, 23.5 to 25.5 mag/D" in J and ~22 to 24 mag/D" in 
K, allow for the detection and study of the NIR structural prop- 
erties and colors of the underlying stellar low-surface bright- 
ness (LSB) component in all sample galaxies. This evolved 



Noeske et al.: NIR observations of Blue Compact Dwarf Galaxies 



27 



LSB host, underlying the star-forming regions is known to ex- 
ist in the majority of BCDs and to contain the bulk of the stellar 
mass in these systems. Consequently, a systematic determina- 
tion of its structural properties (e.g., radial stellar surface den- 
sity distribution, central surface brightness, exponential scale 
length) and of the gravitational potential it forms may prove 
crucial for the understanding of the starburst activity, dynam- 
ics and evolution of BCDs. Other than in optical wavelengths, 
where extended starburst emission hides the LSB host inside 
its inner 2-3 exponential scale lengths, NIR studies allow to 
extend surface photometry of the underlying old stellar back- 
ground to a smaller galactocentric distance, thereby better con- 
strain its overall intensity distribution. Our results can be sum- 
marized as follows: 

1 . Surface brightness profiles (SBPs) of our sample BCDs can 
at large galactocentric radii R* be well approximated by an 
exponential fitting law, in agreement with previous obser- 
vational evidence based on deep optical surface photome- 
try. Also, the exponential scale lengths derived in the opti- 
cal and NIR spectral domain are in mutual agreement, im- 
plying minor optical-NIR color gradients in the LSB com- 
ponent. 

2. The LSB host galaxy of several sample BCDs shows on 
a galactocentric radius comparable to the size of star- 
forming component (within its inner 1-3 exponential scale 
lengths) a conspicuous intensity depression with respect to 
the purely exponential slope, observed for larger radii. This 
type of an inwards flattening exponential SBP, classified 
type V in Binggeli & Cameron (1991), has been frequently 
derived in dwarf irregular and dwarf elliptical galaxies. The 
rare detection of type V profiles in BCDs from previous op- 
tical studies can be attributed to the fact that in those wave- 
lengths extended starburst emission overshines the under- 
lying LSB population within typically its inner ~2 expo- 
nential scale lengths. A possible high frequency of type V 
profiles among dwarf galaxies will not have a notable effect 
on the structural difference between BCDs and dwarf irreg- 
ulars on the parameter space defined by the central surface 
brightness, exponential scale length and absolute magni- 
tude of their LSB host galaxy. However, it may have impor- 
tant implications for our view about BCDs, as it would sig- 
nificantly increase the estimated starburst-to-LSB luminos- 
ity fraction, and therefore the amount of photometric fading 
of these systems, once the starburst activity has terminated. 
This information is crucial for, e.g., establishing or discard- 
ing the hypothesis of faint dwarf spheroidals being the evo- 
lutionary endpoints of BCDs. In the same way, a typeV 
intensity distribution would impose new observational con- 
straints to the derivation of the total stellar mass and its 
radial density distribution within the underlying LSB host 
galaxy of BCDs. 

3. The exact shape and physical origin of typeV SBPs in 
dwarf galaxies are still to be investigated. We find that such 
SBPs can be well approximated by a modified exponen- 
tial fitting formula proposed in Papaderos et al. (1996a). 
Alternatively, a Sersic law can also yield good fits to type V 
profiles, albeit small systematic residuals. However, the 



practical applicability of the Sersic law to the LSB emis- 
sion of BCDs is limited by the strong non-linear coupling 
of its free parameters, and the extreme sensitivity of the 
achieved solutions to, e.g., small uncertainties in the sky 
subtraction and SBP derivation. 

4. For the majority of the LSB host galaxies in our sample 
we derive optical-NIR colors indicative of an evolved stel- 
lar population of subsolar metallicity. The metal-deficient 
BCDs Tol 65 and Tol 1214-277 show blue colors (B - 
J <1.3 mag) in their underlying stellar LSB component. 

5. NIR images, in combination with optical data reveal for 
some BCDs signatures of appreciable and non-uniform 
dust absorption on a spatial scale as large as ~1 kpc. The 
NIR images, being less affected by dust extinction and 
widespread ionized gas emission than optical data, allow 
for the detection of a variety of morphological features 
within the star-forming component, as e.g. coherent as- 
semblies of compact stellar clusters, that might be rem- 
nants from previous episodes of collective star formation. 
An elaborate multi wavelength investigation of such sources 
holds the promise of a better understanding of the history 
and spatial progression of star-forming activity in BCDs. 

Acknowledgements. Research by K.G.N, has been supported by 
the Deutsche Forschungsgemeinschaft (DFG) Grants FR325/50-1 
and FR325/50-2. LMC acknowledges support from the European 
Community Marie Curie Grant HPMF~CT-2000-00774. We thank 
U. Fritze - v. Alvensleben, R Anders, J. Bicker and J. Schulz for 
kindly providing the GALEV models. This research has made use 
of the NASA/IR\C Extragalactic Database (NED) which is oper- 
ated by the Jet Propulsion Laboratory, CALTECH, under contract 
with the National Aeronautic and Space Administration. This pub- 
lication makes use of data products from the Two Micron All Sky 
Survey, which is a joint project of the University of Massachusetts 
and the Infrared Processing and Analysis Center/California Institute 
of Technology, funded by the National Aeronautics and Space 
Administration and the National Science Foundation. 

References 

Alton, R B., Draper, R W., Gledhill, T. M., et al. 1994, MNRAS, 270, 
238 

Andredakis, Y. C, Peletier, R. F & Balcells, M. 1995, MNRAS, 275, 
874 

Andreon, S. 2002, A&A, 382, 821 

Babul, A. & Rees, M. J. 1992, MNRAS, 255, 346 

Beck, S. C, Kelly, D. M. & Lacy, J. H. 1997, AJ, 1 14, 585 

Bergvall, N. & Ostlin, G. 2002, A&A, 390, 891 

Bertelli, G., Bressan, A., Chiosi, C, Fagotto, F. & Nasi, E. 1994, 

A&AS, 106, 275 
Binggeh, B. & Cameron, L. M. 1991, Ac&A, 252, 27 
Binggeh, B. & Jerjen, H. 1998, A&A, 333, 17 
Cairos, L. M. 2000, PhD Thesis, Instituto de Astrofisica de Canarias 
Cairos, L. M., Caon, N., Garcia-Lorenzo, B., Vflchez, J. M. & Muiioz- 

Tuiion, C. 2002a, ApJ, 577, 164 
Cairos, L. M., Caon, N., Papaderos, P., et al. 2003, ApJ, in press (C03) 
Cairos, L. M., Caon, N., Vflchez, J. M., Gonzalez-Perez, J. & Muiioz- 

Tunon, C. 2001b, ApJS, 136, 393 
Cairos, L. M., Vflchez, J. M., Gonzalez Perez, J., Iglesias-Paramo, J. 

& Caon, N. 2001a, ApJS, 133, 321 
Campbell, A. W. & Terlevich, R. 1984, MNRAS, 21 1, 15 



28 



Noeske et al.: NIR observations of Blue Compact Dwarf Galaxies 



Cannon, J. M., Skillman, E. D., Gamett, D. R. & Dufour, R. J. 2002, 
ApJ, 565, 931 

Caon, N., Capaccioli, M. & D'Onofrio, M. 1993, MNRAS, 265, 1013 

Cardelli, J. A., Clayton, G. C. & Mathis, J. S. 1989, ApJ, 345, 245 

Carpenter, J. M. 2001, AJ, 121, 2851 

Cellone, S. A. & Buzzoni, A. 2001, A&A, 369, 742 

Cellone, S. A., Forte, J. C. & Geisler, D. 1994, ApJS, 93, 397 

Conti, R S. & Vacca, W. D. 1994, ApJ, 423, L97 

Corbin, M. R., Korista, K. T. & Vacca, W. D. 1993, AJ, 105, 1313 

Cutri, R. M., et al. 2000, The 2MASS Explanatory Supplement 

Davidge, T.J. 2002, PASR 114, 1206 

Davies, J. I. & Phillipps, S. 1988, MNRAS, 233, 553 

Dekel, A. & Silk, J. 1986, ApJ, 303, 39 

de Vaucouleurs, G., de Vaucouleurs, A., Corwin, H. G., et al. 1991, 

Third Reference Catalogue of bright Galaxies (Springer Verlag) 
DouWier, V., Caulet, A. & Comte, G. 1999, A&AS, 138, 213 
Doublier, V., Caulet, A. & Comte, G. 2001, Ac&A, 367, 33 
Doublier, V., Comte, G., Petrosian, A., Surace, C. & Turatto, M. 1997, 
A&AS, 124, 405 

Due, R-A., Papaderos, R, Balkowski, C, et al. 1999, A&AS, 136, 539 
Elias, J. H., Frogel, J. A. & Humphreys, R. M. 1985, ApJS, 57, 91 
Ferrara, A. & Tolstoy, E. 2000, MNRAS, 313, 291 
Fioc, M. & Rocca-Volmerange, B. 1997, A&A, 326, 950 
Freeman, K. C. 1970, ApJ, 160, 811 

Fricke, K. J., Izotov, Y. I., Papaderos, P., Guseva, N. G. & Thuan, T. 

X. 2001, AJ, 121, 169 (FOl) 
Fritz, T. 2000, PhD Thesis, University of Bonn 
Fujita, A., Mac Low, M.-M., Meiksin, A. & Ferrara, A. 2001, 

American Astronomical Society Meeting, 198, 5703 
Gallagher, J. S. III. & Hunter, D. A. 1989, AJ, 98, 806 
Garcia, A. M. 1993, A&AS, 100, 47 

Gil de Paz, A., Zamorano, J., Gallego, J. & Domfnguez, F. de B. 

2000a, A&AS, 145, 377 
Gil de Paz, A., Zamorano, J. & Gallego, J. 2000b, A&A, 361, 465 
Gonzalez-Riestra, R., Rego, M. & Zamorano, J. 1987, A&A, 202, 72 
Graham, A., Lauer, T. R., Colless, M. & Postman, M. 1996, ApJ, 465, 

534 

Guseva, N. G., Izotov, Y. I., Papaderos, P, et al. 2001, A&A, 378, 756 
Guseva, N.G., Izotov, Y.I. & Thuan, TX. 2000, ApJ, 531, 776 
Guzman, R., Jangren, A., Koo, D. C, Bershady, M. A. & Simard, L. 

1998, ApJ, 495, L13 
Hensler, G., Dickow, R. & Junkes, N. 1997, RMxAA Conf. Series, 6, 

90 

Heydari-Malayeri, M., Melnick, J. & Martin, J.-M. 1990, A&A, 234, 
99 

Hidalgo-Gamez, A. M., Masegosa, J. & Olofsson, K. 2001, A&A, 
367, 388 

Hidalgo-Gamez, A. M., Masegosa, J. & Olofsson, K. 2002, RMxAA 

Conf. Series, 12, 242 
Hirashita, H., Tamura, N. & Takeuchi, T. T. 2000, lAU Symposium, 

201, 64 

Hunt, L.K., Thuan, TX. & Izotov, Y.I. 2003, ApJ, 588, 281 

Hunter, D. & Hoffman, L. 1999, AJ, 117, 2789 

Izotov Y. I., Chaffee, F H. & Green, R. F 2001, ApJ, 562, 727 

Izotov, Y. I. & Thuan, T. X. 1999, ApJ, 51 1, 639 

Izotov, Y. I. & Thuan, T. X. 2002, ApJ, 567, 875 

Izotov, Y. I., Thuan, T. X. & Lipovetsky, V. A. 1997, ApJS, 108, 1 

James, P A. 1994, MNRAS, 269, 176 

Jarrett, T. H., Chester, T, Cutri, R., Schneider, S. & Huchra, J. 2000, 
AJ, 119,2498 

Johnson, K. E., Leitherer, C, Vacca, W. D. & Conti, P S. 2000, AJ, 
120, 1273 

Karachentsev, I. D., Sharina, M. E., Dolphin, A. E., et al. 2002, A&A, 
385, 21 



Kobulnicky, H. A., Kennicutt, R. C. & Pizagno, J. L. 1999, ApJ, 514, 
544 

Kriiger, H. 1992, PhD Thesis, University of Gottingen 
Kriiger, H., Fritze-v. Alvensleben, U. & Loose, H.-H. 1995, A&A, 
303,41 

Kunth, D., Maurogordato, S. & Vigroux, L. 1988, A&A, 204, 10 

Kunth, D. & Sargent, W. L. W. 1981, A&A, 101, L5 

Kunth, D. & Sargent, W. L. W. 1983, ApJ, 273, 81 

Lin, D. N. C. & Faber, S. M. 1983, ApJ, 266, L21 

Loose, H. H. & Thuan, T. X. 1986, Star Forming Dwarf Galaxies and 

Related Objects, 73 (LT86) 
Lowenthal, J. D., Koo, D. C, Guzman, R., et al. 1997, ApJ, 481, 673 
MacArthur, L. A., Courteau, S. & Holtzman, J. A. 2002, ApJ, in press 

lastro-ph/0208404i 
Makarova, L. N., Karachentsev, I. D., Grebel, E. K. & Barsunova, O. 

Y. 2002, A&A, 384, 72 
Makarova, L., Karachentsev, I., Takalo, L. O., Heinaemaeki, P. & 

Valtonen, M. 1998, A&AS, 128, 459 
Marlowe, A. T, Heckman, T. M., Wyse, R. F. G. & Schommer, R. 

1995, ApJ, 438, 563 
Marlowe, A. T, Meurer, G. R. & Heckman, T. M. 1999, ApJ, 522, 183 
Marlowe, A. T, Meurer, G. R., Heckman, T. M. & Schommer, R. 

1997, ApJS, 112, 285 
Masegosa, J., Moles, M. & Campos-Aguilar, A. 1994, ApJ, 420, 576 
Mas-Hesse, J. M. & Kunth, D. 1999, A&A, 349, 765 
Mendez, D.I. & Esteban, C. 1999, AJ, 118, 2723 
Mendez, D.I. & Esteban, C. 2000, A&A, 359, 493 
Meurer, G. R., Staveley-Smith, L. & Killeen, N. E. B. 1998, MNRAS, 

300, 705 

Mirabel, I. F & Sanders, D. B. 1988, ApJ, 335, 104 
Noeske, K. G. 1999, Diploma Thesis, University of Gottingen 
Noeske, K. G., Guseva, N. G., Fricke, K. J., et al. 2000, A&A, 361, 31 
Noeske, K. G., Iglesias-Paramo, J., Vflchez, J. M., Papaderos, P. & 

Fricke, K. J. 2001, A&A, 371, 806 
Noguchi, M. 2001, ApJ, 555, 289 
Ott, J., Walter, F, Brinks, E., et al. 2001, AJ, 122, 3070 
Pagel, B. E. J., Simonson, E. A., Terlevich, R. J. & Edmunds, M. G. 

1992, MNRAS, 255, 325 
Papaderos, P., et al. 2003, in prep. 

Papaderos, P. 1998, PhD Thesis, University of Gottingen 
Papaderos, P. & Fricke, K. J 1998, Highlights in X-ray Astronomy, B. 

Aschenbach et al. (eds.) 
Papaderos, P., Fricke, K. J., Thuan, T. X., Izotov, Y. I. & Nicklas, H. 

1999, A&A, 352, L57 (P99) 
Papaderos, P, Izotov, Y.I., Thuan, T.X., et al. 2002, A&A, 393, 461 
Papaderos, P, Loose, H.-H., Fricke, K. J. & Thuan, T. X. 1996b, A&A, 

314,59 (P96b) 

Papaderos, P, Loose, H.-H., Thuan, TX. & Fricke, K.J. 1996a, 

A&AS, 120, 207 (P96a) 
Pastoriza, M. G. & Dottori, H. A. 1981, ApJ, 244, 27 
Patterson, R. J. & Thuan, T. X. 1996, ApJS, 107, 103 
Peimbert, M. & Torres-Peimbert, S. 1974, ApJ, 193, 327 
Persson, S. E., Murphy, D. C, Krzeminski, W., Roth, M. & Rieke, M. 

L 1998, AJ, 116, 2475 
Pustilnik, S. A., Kniazev, A. Y, Lipovetsky, V. A. & Ugryumov, A. V. 

2001, A&A, 373, 24 
Rieschick, A. & Hensler, G. 2001, ApSS SuppL, 277, 111 
Sage, L. J., Salzer, J. J., Loose, H.-H. & Henkel, C. 1992, A&A, 265, 

19 

Salzer, J. J., MacAlpine, G. M. & Boroson, T. A. 1989, ApJS, 70, 479 
Salzer, J. J. & Norton, S. A. 1999, ASP Conf. Ser. 170: The Low 

Surface Brightness Universe, 253 
Schaerer, D., Contini, T. & Kunth, D. 1999, A&A, 341, 399 



Noeske et al.: NIR observations of Blue Compact Dwarf Galaxies 



29 



Schlegel, D. J., Finkbeiner, D. P. & Davis, M. 1998, ApJ, 500, 525 
Schulte-Ladbeck, R. E., Hopp, U., Greggio, L. & Crone, M. M. 2000, 
AJ, 120, 1713 

Schulte-Ladbeck, R. E., Hopp, U., Greggio, L., Crone, M. M. & 

Drozdovsky, I. O. 2001, AJ, 121, 3007 
Schulz, J., Fritze - v. Alvensleben, U., MoUer, C. S. & Fricke, K. J. 

2002, A&A, 392, 1 
Sersic, J.-L. 1968, Atlas de Galaxias Australes, Observatorio 

Astronomico de Cordoba 
Silk, J., Wyse, R. F G. & Shields, G. A. 1987, ApJ, 322, L59 
Simpson, C. E. & Gottesman, S. T. 2000, AJ, 120, 2975 
Smoker, J. V., Davies, R. D., Axon, D. J. & Hummel, E. 2000, A&A, 

361, 19 

Stasinska, G., Comte, G. & Vigroux, L. 1986, A&A, 154, 352 

Tammann, G. A. & Sandage, A. 1985, ApJ, 294, 81 

Taylor, C. L., Brinks, E., Grashuis, R. M. & Skillman, E. D. 1995, 

ApJS, 99, 427 
Telles, E. & Terlevich, R. 1995, MNRAS, 275, 1 
Telles, E., Melnick, J. & Terlevich, R. 1997, MNRAS, 288, 78 
Thuan, T. X. 1991, Massive Stars in Starbursts, C. Leitherer et al. 

(eds.) 

Thuan, T X. 1985, ApJ, 299, 881 

Thuan, T. X., Izotov, Y. 1. & Foltz, C. B. 1999, ApJ, 525, 105 
TuUy, R. B. 1988, Nearby Galaxies Catalog (Cambridge University 
Press) 

Vacca W.D. & Conti, P.S. 1992, ApJ, 401, 543 
Vallenari, A. & Bomans, D. J. 1996, A&A, 313, 713 
van den Bergh, S. 2000, AJ, 119, 609 
van Zee, L. 2000, AJ, 1 19, 2757 

van Zee, L., Skillman, E.D. & Salzer, J.J. 1998, AJ, 116, 1186 
van Zee, L., Salzer, J.J. & Skillman, E.D. 2001, AJ, 122, 121 
Vanzi, L., Hunt, L. K. & Thuan, T. X. 2002, A&A, 390, 481 
Vanzi, L., Hunt, L. K., Thuan, T. X. & Izotov, Y. I. 2000, A&A, 363, 
493 

Vanzi, L., Rieke, G. H., Martin, C. L. & Shields, J. C. 1996, ApJ, 466, 
150 

Vennik, J., Hopp, U. & Popescu, C. C. 2000, A&AS, 142, 399 
Vigroux, L., Souviron, J., Lachieze-Rey, M. & Vader, J. P. 1988, 
A&AS, 73, 1 

Wainscoat, R. J. & Cowie, L. L., 1992, AJ, 103, 332 
Walter, F, Brinks, E., Duric, N. & Klein, U. 1997, AJ, 113, 2031 
Yasuda, N., Fukugita, M. & Okamura, S. 1997, ApJS, 108, 417 
Young, C. K. & Currie, M. J. 1994, MNRAS, 268, Lll 
Young, C. K. & Currie, M. J. 1995, MNRAS, 273, 1141 



